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Prediction  of  Aerodynamic  Drag 


1  INTRODUCTION 

The  state-of-the-art  in  aircraft  drag  prediction  was  defined  as  of  1973 

20 

by  the  AGARD  conference  entitled  "Aerodynamic  Drag."  This  chapter  draws 
heavily  on  that  work  and  attempts  to  update,  to  the  extent  possible,  those 
portions  relevant  to  aircraft  drag  prediction  at  subsonic,  transonic  and 
supersonic  speeds.  Special  drag  prediction  problems  peculiar  to  short  take¬ 
off  and  landing  (STOL)  designs,  such  as  jet  flap  thrust  recovery,  and  energy 
efficient  aircraft  concepts,  for  example  laminar  flow  control,  are  beyond  the 
scope  of  this  chapter;  as  are  drag  due  to  speed  brakes  and  fighter  aircraft 
weapons  carriage.  Certainly,  a  sound  physical  understanding  of  drag  will  be 
required  before  drag  prediction  becomes  a  science. 

The  assessments,  projections  and  conclusions  in  this  chapter  agree,  in 
general,  with  those  of  Wood161  who  briefly  addressed  four  relevant  questions. 

"1.  How  well  do  ground  based  estimates  for  drag  polars  and  engine  charac¬ 
teristics  correlate  to  flight  test  results? 

2.  In  what  areas  do  performance  prediction  techniques  work  best  or 
worst,  and  why? 

3.  What  are  the  differences  in  the  way  various  manufacturers  (airframe 
or  engine)  predict  drag  polars  and  engine  characteristics? 

4.  Are  there  portions  of  various  performance  prediction  techniques  that 
could  or  should  be  combined  to  form  a  better  or  best  method?" 
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Wood's  answer  to  question  one  is,  very  well  for  up-and-away  flight; 
however,  his  survey  of  twelve  commercial  transport  aircraft  showed  six 
predictions  were  low  (as  much  as  22%),  four  were  high  (up  to  10%),  and  two 
were  exact.  The  second  question  received  a  similar  answer:  predictions  are 
best  for  up-and-away  flight  in  the  clean  aircraft  configuration.  An  entire 
symposium  would  be  needed  to  address  questions  three  and  four  unless  the 
simple  answer  that  no  single  methodology  is  best  for  every  situation  is 
accepted.  Wood  also  observed  that,  "the  quality  of  the  estimate  is  more  a 
function  of  the  time  and  care  taken  to  include  all  the  details  and  higher- 
order  terms  than  it  is  of  the  particular  equations  used." 

1 I  CQ 

Williams'  synthesis  of  responses  to  the  AGARD-FMP  questionnaire  on 
prediction  techniques  and  flight  correlation  is  also  worthy  of  detailed  study. 
The  five  questions  asked  of  experts  in  six  NATO  countries  were: 

"1.  What  are  the  advantages/disadvantages  of  different  prediction  tech¬ 
niques? 

2.  What  portions  of  the  flight  regime  cannot/should  not  be  addressed  by 
ground-based  techniques? 

3.  Are  there  areas  where  analytical  prediction  can  be  better  than  wind 
tunnel  and/or  simulation  results;  or  vice-versa? 

4.  Are  there  methods  of  reducing  differences  between  prediction  and 
flight  test  results? 

5.  Are  there  any  new  prediction  techniques  that  should  be  emphasized?" 

Williams'  section  on  background  to  prediction/design  needs  and  discussion 
of  increasing  technical  demands  is  also  especially  relevant  to  the  understand¬ 
ing  of  drag  prediction. 
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The  initial  sections  of  this  chapter  describe  drag  prediction  methods 
typically  used  to  define  the  most  promising  configuration  for  further,  de¬ 
tailed  analysis  and  wind  tunnel  testing.  These  methods  are  mainly  empirical 
in  origin,  although  numerical  aerodynamics  is  presently  providing  timely  data 
for  this  phase  of  design  in  some  organizations.  The  latter  sections  of  the 
report  describe  drag  prediction  methods  based  on  wind  tunnel  results  from 
models  specifically  constructed  for  the  particular  design  project  and  guided 
by  aerodynamic  theory,  known  as  "wind  tunnel  to  flight  correlations." 
Rooney11^*^  and  Craig3^  have  recently  published  detailed  accounts  of  the 
meticulous  testing  performed  to  correlate  the  wind  tunnel  and  flight  measured 
aerodynamic  drag  of  the  Tomahawk  Cruise  Missile.  The  size  of  this  missile 
permitted  wind  tunnel  testing  of  the  full-scale  vehicle  --  yet  discrepancies 
remained.  These  correlations  and  many  others  for  all  types  of  aircraft  and 
missiles  form  the  basis  of  evaluation  for  a  new  design. 

Despite  many  recent  advances,  it  is  generally  conceded  that  accurate  drag 
predictions,  based  entirely  on  the  solution  of  the  equations  of  motion,  or 
some  reduced  form  thereof;  are,  at  present,  limited  to  the  prediction  of  drag 
due  to  lift  in  the  linear  range  and  supersonic  wave  drag  for  a  limited  class 
of  slender  configurations.  All  other  drag  predictions  ultimately  depend  on 
empirical  correlations.  Progress  in  computational  fluid  dynamics  or,  more 
generally,  numerical  aerodynamics,  has  been  tremendous  in  the  last  decade  and 
will  continue  into  the  future  due  to  the  introduction  of  new  computers, 
faster,  more  accurate  solution  algorithms,  improved  resolution  of  grids  and 
new  turbulence  models.  However,  except  for  the  isolated  cases  of  drag-due- 
to-lift  at  small  angle  of  attack  and  supersonic  wave  drag  for  smooth  slender 
bodies,  drag  prediction  is  beyond  the  capability  of  current  numerical  aero¬ 
dynamic  methods. 


3 


The  aerodynamics  of  the  aerodynamic  reference  model usually  a  sub¬ 
scale  wind  tunnel  wing-body  model  with  flow-through  nacelles  and  reference 
nozzles  are  considered  initially.  This  is  one  of  many  models  tested  during  an 
aircraft  development  program.  Fig.  1,  and  the  basic  model  for  drag  prediction. 
The  uncorrected  airplane  lift,  drag,  and  pitching  moment  characteristics  are 
derived  from  tests  of  this  model.  Drag  corrections  to  account  for  the  differ¬ 
ences  between  the  data  from  this  model  and  the  actual  flight  vehicle  are  dis¬ 
cussed  in  the  latter  sections  of  this  chapter.  The  methods  of  properly 
combining  these  drag  data,  including  thrust  effects,  are  the  subject  of  this 
book. 

1.1  Drag 

The  resultant  aerodynamic  force  caused  by  a  flight  vehicle's  motion  with 
respect  to  the  atmosphere  is  the  summation  of  the  normal,  pressure  forces,  and 
the  tangential  skin  friction  forces  acting  on  the  vehicle's  surface.  This 
resultant  force  is  conventionally  resolved  into  lift  and  drag  components  in 
the  vehicle's  plane  of  symmetry.  The  lift  force  is  the  aerodynamic  reaction 
perpendicular  to  the  free-stream  velocity  direction  which  is  the  level 
equilibrium  flight  path  direction.  Lift  is  not  directly  a  subject  of  this 
chapter;  however,  drag,  the  component  of  the  total  force  that  opposes  motion 
in  the  equilibrium  flight  path  direction,  approximately  follows  a  parabolic 
variation  with  lift  for  heavier-than-air  flight  as  shown  in  Figure  2.  Thus, 
the  relative  state  of  motion  of  the  vehicle  (equilibrium,  or  accelerated 
flight)  is  determined  by  the  lift  force,  which  is  an  appropriate  multiple  of 
the  vehicle's  weight,  and  the  throttle-dependent  thrust  force. 

Aircraft  normally  spend  ninety  percent  of  their  flight  time  on  nominally 
straight,  unaccelerated  flight  paths  where  all  the  forces  are  in  static 
equilibrium.  This  is  the  cruise  condition  that  is  usually  considered  the 
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standard  condition  for  the  design  of  an  airplane.  This  chapter  is  mainly 
concerned  with  cruise  drag  prediction  (Fig.  2)  because  of  the  importance  of 
this  flight  condition.  Drag  prediction  for  accelerated  flight  (take-off, 
landing,  maneuver)  will  be  treated  as  an  important,  although  ancilliary, 
issue. 

The  parabolic  variation  of  drag  with  lift,  or  angle  of  attack,  appropri¬ 
ately  additive  to  the  zero  lift  drag  is  shown  in  Fig.  2.  Drag  is  produced  by 
the  tangential  (skin  friction)  and  normal  (pressure)  forces  acting  on  the 
vehicle's  surface  due  to  relative  fluid  motion  from  a  basic  fluid  mechanics 
standpoint.  Total  drag  has  been  dissected  into  a  multitude  of  various  compo¬ 
nents,  depending  on  the  phase  of  the  design  process,  the  aircraft  mission  and 
configuration,  the  experimental  data  available  and  the  persuasion  of  the  drag 
analyst.  A  properly  defined,  consistent  reference  condition  for  the  entire 
aircraft  and  propulsion  system  is  crucial  to  the  correct  final  comparison  of 
theoretical,  empirical,  ground  and  flight  test  data  in  the  latter  design 
stages  for  performance  guarantees.  The  operating  reference  condition  is 
straight-forward  for  commercial  transports  with  a  single,  known  cruise  design 
point;  however,  military  fighter  aircraft  with  multiple  design  points  pose 
severe  additional  problems.  It  is  often  difficult  to  define  an  operating 
reference  condition  in  this  case.  Additionally,  the  variation  in  installed 
flight  test  thrust  permitted  in  military  engines  of  the  same  family  can  be  an 
order  of  magnitude  greater  than  that  permitted  in  engines  slated  for 
commercial  applications. 

1.2  Level  of  Drag  Prediction  Detail 

Three  somewhat  distinct  levels  of  drag  prediction  sophistication  and 

reliability  are  usually  described  by  authors  on  this  subject.20,  106,  110, 
125167 

9  In  reality,  drag  predictions  are  updated  continuously  throughout  the 
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entire  design  cycle  of  an  aircraft  as  data  from  numerical  methods  and  wind 
tunnel  tests  become  available.  However,  decisions  that  affect  the  final 
design  are  made  based  on  the  best  available  information  at  many  stages  of  the 
design  process. 

1.2.1  The  Preliminary  Design  Phase. 

This  is  the  beginning  of  the  design  process.  It  has  also  been  called  the 

feasibility^®  and  conceptual*®'7  design  phase.  Primarily  empirical  methods  are 

used  to  assess  the  relative  merits  of  many  design  concepts  against  the  mission 

specifications  generated  by  an  apparent  market  opportunity  or  military 

statement  of  operational  need.  Drag  prediction  error,  as  measured  with 

respect  to  future  flight  test  data,  is  highest  due  to  method  error  and 

13 

geometric  uncertainty  in  the  configuration  definition  ,  Fig.  3.  Relative 

accuracy  in  the  methodology  is  necessary  at  this  stage  in  order  to  select  the 

most  efficient  configuration  to  meet  the  design  objectives.  Absolute  accuracy 

is  desirable  in  order  to  give  equal  consideration  to  all  new  concepts  and 

areas  of  the  design  space  or  mission  profile.  Promising  new  concepts  could  be 

excluded  from  further  consideration  at  this  stage  if  drag  methodology 

incorrectly  predicts  that  the  mission  cannot  be  achieved.  An  overall  drag 

target  or  maximum  drag  level  is  determined  during  this  design  phase.  Several 

compilations  of  drag  prediction  methodology  are  available^’  *®*’  *®7,  *2^’ 
135 

for  rapid  performance  assessment.  Each  has  its  own  peculiarities  and 
limitations.  Additionally,  each  airframe  manufacturer  has  compiled  large  drag 
handbooks  that  are  highly  valued  and  extremely  proprietary. 

1.2.2  The  Detailed  Design  Phase. 

Drag  validation  becomes  a  more  specialized  and  detailed  continuous  process. 

This  phase  has  also  been  called  the  development^®,  preliminary*^,  and  project 
109 

definition  phase.  Estimates  and  assumptions  are  gradually  replaced  by 
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predictions  supported  by  intensive  wind  tunnel  test  data  on  the  determinate 
aircraft  design,  as  derived  from  the  clean  wing  reference  configuration  and 
aerodynamic  reference  configuration  models,  Fig.  1.  Drag  sub-targets  for  each 
component  of  the  configuration,  consistent  with  the  overall  target,  are 
assessed  and  synthesized.  In  some  development  programs,  prototype  flight  test 
data  provides  the  basis  for  diagnostics,  reassessments,  and  analysis  if  drag 
improvements  are  found  to  be  necessary. 

1.2.3  The  Final  Design  Phase. 

Production  aircraft  performance  guarantees  with  error  bounds  are  made. 

This  is  also  the  pre-production^®  or  wind  tunnel  to  flight  correlation1®4  phase. 
A  detailed  analysis  and  interpretation  of  the  prototype  test  performance  is 
conducted.  The  drag  predictive  model  is  calibrated  by  reference  to  the 
prototype  data.  An  excellent  description  of  this  process  and  its  difficulties 
has  been  published  by  Rooney119  and  Craig34  for  the  Tomahawk  missile. 

Theory,  empirical  methods,  and  ground  and  flight  test  measurements  are 
continuously  intertwined  throughout  the  design,  development  and  production 
phases  of  an  aircraft's  life.  Each  drag  prediction  method  contains  its  unique 
strengths  and  weaknesses. 

Strictly  empirical  methods,  conventionally  based  on  flight  test  data, 
include  the  real  aircraft  effects  of  flow  separations,  Reynolds  number,  gaps, 
steps,  protuberances,  rigging,  etc.,  but  are  not  readily  generalized.  Prob¬ 
lems  have  even  arisen  when  these  methods  are  applied  to  aircraft  of  the  same 
family  from  which  the  data  base  was  generated. 

Current  theoretical  methods  that  utilize  large  computers  are  capable  of 
immediate  generalization  to  almost  any  configuration,  but  are  limited  in 
accuracy  because  of  the  lack  of  real  aircraft  effects.  Prediction  of  the 
drag-due-to-lift  is  usually  quite  accurate  up  to  drag  levels  where  flow 
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separations  become  significant.  Prediction  of  the  zero-lift-drag  remains 
empi rical . 

Progressive  wind  tunnel  development  programs  contain  a  greater  degree  of 

reality  than  theoretical  methods,  lack  some  real  aircraft  effects,  and  are 

costly  and  time  consuming.  Correct  interpretation  of  the  data  remains  an  art, 

particularly  when  wall  effects,  Reynolds  number  scaling  and  propulsion  system 

model  data  are  included.  Systematic  total  drag  measurement  errors  greater 

than  one  percent  cannot  be  tolerated  if  drag  prediction  methods  for  the 

19 

smaller  drag  components  are  to  be  validated. 

A  typical  breakdown  of  flight  test  drag  data  is  shown  in  Fig.  2.  The 
industry  standard  for  drag  measurement  (±  3%  accuracy)  is  seldom  achieved, 

IIO 

even  under  carefully  controlled  conditions.  From  these  data,  drag 

correlations  for  each  drag  component  are  made  and  drag  design  charts  are 

constructed  as  a  basis  for  the  preliminary  design  empirical  method.  The 

preliminary  design  study  is  considered  as  the  reference  point  for  judgment  on 

program  success.  The  initial  decisions  on  wing  area,  aspect  ratio,  thickness 

ratio  and  matching  engine  cycle  and  thrust  requirements  are  based  on  early 

lift-to-drag  (L/D)  estimates.  A  lower  bound  for  L/D  must  be  estimated  in 

order  to  assure  that  the  propulsion  system  will  provide  the  power  necessary  at 

the  airplane  critical  flight  condition.  Misjudgements  at  this  design  stage 

will  result  in  an  inefficient  engine-airplane  that  is  noncompetitive. 

The  designer  encounters  errors  from  two  sources  at  this  stage;  errors  due 

to  lack  of  configuration  geometric  definition,  and  errors  due  to  inaccurate 

13 

methods.  These  errors  have  been  estimated  by  Bowes  for  long-range  subsonic 
transports  and  are  summarized  in  Fig.  3.  He  also  suggests  an  aerodynamicist's 
rating  scale,  where  M^  is  the  drag  divergence  Mach  Number.  The  need  for  an 
accurate  force  and  moment  model  as  early  as  possible  in  the  design  cycle  is 
clearly  evident,  from  the  following  tables. 
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L/D  Level  Achieved 

mdd 

Rating 

±  3% 

±  .002 

Amazing 

±  5% 

±  .004 

Very  Good 

±  7% 

±  .006 

Average 

±  10% 

±  .010 

Below  Average 

Most  commercial  aircraft  flying  today  fit  within  this  table;  however,  military 
programs  have  been  far  less  successful  in  achieving  these  accuracy  levels. 

This  may  be  due,  in  part,  to  the  less  precise  flight  test  data  recorded  from 
the  older  military  aircraft  as  shown  in  the  following  table. 

TABLE  1  DRAG  SCATTER 


Aircraft 

Flight 

Test  Data  Scatter 

Ref 

B-747 

±  1% 

(drag) 

13 

C-5A 

±  3.5%  (CDp)  ±  2%  (drag) 

86 

C- 141 

±  3.3%  (drag) 

86 

B-707 

±  5% 

(drag) 

13 

DC-10-30 

±  3.5% 

(range  factor) 

44 

F-14 

±  1% 

118 

F— 15 

±  20% 

(drag) 

173 

YF-16 

±  5% 

(drag) 

18 

F— 16 

±  12% 

^Dmin^ 

152 

F- 18 

±  20% 

(drag) 

Tomahawk 

±1.5% 

(drag) 

34 

Tornado 

±  3% 

(drag)  ±  1%  (CDq) 

71 

Alpha-Jet 

±  5% 

(drag) 

124 
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The  table  should  not  be  read  such  that  the  data  scatter  is  constant  for  a 
given  aircraft  throughout  the  flight  envelope.  The  scatter  stated  is 

3 

representative  of  the  cruise  condition.  Fig.  4,  from  Arnaiz's  study  of  the 
XB-70  flight  test  data,  clearly  demonstrates  this  point  by  showing  the  actual 
number  of  occurrences  of  scatter  and  the  estimated  scatter  plotted  against 
Mach  number.  The  drag  uncertainty  estimates  were  obtained  by  combining  many 
individual  errors  in  weight,  Mach  number,  static  pressure,  angle  of  attack, 
thrust,  acceleration,  etc.,  in  a  root-sum-square  calculation. 

Performance  estimates  for  military  aircraft  are  more  accurate  than  Table 
1  indicates,  since  the  major  problems  occur  when  one  attempts  to  isolate  the 
thrust  and  drag  components  from  the  measured  performance  data.  Current 
steady-state,  level  flight  test  data  indicate  that  a  repeatability  of  ±  1%  can 
be  achieved  if  meticulous  care  is  exercised  during  the  engine  calibration  and 
flight  test  processes.  This  level  of  accuracy  must  be  attained159  in  order  to 
obtain  useful  correlations  for  the  smaller  drag  components  shown  in  Fig.  2. 
These  component  correlations,  necessary  to  develop  accurate  prediction 
methods,  are  based  on  convention,  tradition,  and  aerodynamic  theory. 

Flight  test  data  are,  ho wever,  not  presently  capable  of  resolving  drag 
into  components.  The  accuracy  of  each  component  cannot  be  checked  since  only 
total  drag,  deduced  from  performance  measurements,  is  available.  Additional 
instrumentation  that  could  resolve  drag  into  components  is  superfluous  to 
aircraft  performance  guarantees  and  is  usually  not  included  on  the  flight  test 
aircraft  because  of  the  additional  cost  to  the  development  program.  This  is 
particularly  true  of  military  aircraft  development  programs  where  budget  cuts 
are  traditional  and  there  is  little  interest  in  developing  drag  prediction 
methodology  for  the  next  aircraft  program. 
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The  relative  size  of  each  drag  component  varies  with  aircraft  type  and 

flight  condition.  Representative  examples  of  the  cruise  drag  and  take-off 

drag  breakdowns  for  four  different  aircraft  types  are  shown  in  Fig.  5  from 
19 

Butler  .  The  drag  breakdown  also  varies  with  aircraft  range  and  purpose 

(military  or  civilian  transport)  within  each  aircraft  type. 

The  conventional  approach  to  drag  estimation  is  to  estimate  and  sum  the 

zero  lift  drag  of  each  major  physical  component  of  the  aircraft  with  allowances 

for  interference  effects  and  other  small  contributions  and  then  add  the  drag 

due  to  lift.  Their  summation  is  intrinsically  inaccurate  since  aircraft 

components  are  mounted  at  incidence  to  each  other  so  each  is  not  at  zero  lift 

19 

simultaneously.  Butler  suggests  that  this  classification  of  aircraft  drag 
into  components  independent  of,  and  dependent  on,  lift  should  be  rejected  and 
drag  could  be  estimated  by  compounding  elements  arising  from  different  basic 
causes  associated  with  fluid  dynamics  as  shown  in  Fig.  6. 

Currently,  each  manufacturer  has  a  methodology,  or  drag  handbook,  by 
which  a  thorough  drag  estimate  is  made.  There  is  no  single,  industry-wide 
drag  handbook.  Manufacturers  may  even  use  several  drag  accounting  breakdowns 
(depending  on  the  customer)  for  the  same  airplane.  There  exist  several  basic 
sources  for  fundamental  drag  level  prediction;  however,  total  drag  compari¬ 
sons,  between  prediction  and  test,  require  accurate  flight  data  and  expert 

I  Cl 

analysis  and  are  highly  proprietary  to  each  manufacturer. 

2  Subsonic  Drag 

The  aircraft  drag  at  subsonic  speeds,  here  taken  to  be  less  than  the  drag 
divergence  Mach  number,  is  traditionally  decomposed  into  lift-induced  drag  and 
minimum  drag.  Minimum  drag  is  further  divided  into  friction  drag,  profile 
drag,  and  interference  drag. 
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The  accuracy  of  the  drag  divergence  Mach  number,  prediction  varies 

119 

widely  with  configuration.  Rooney's  detailed  assessment  of  the  Tomahawk 
flight  test  data  indicates  that  cannot  be  accurately  predicted  while 

Henne,  Dahlin  and  Peavey  conclude  that  a  substantial  capability  exists  for 
the  determination  of  drag  divergence  Mach  number  for  transport  type  geome¬ 
tries.  The  determination  of  drag  levels  and  increments  is  rated  marginal. 

94 

McGeer  and  Shevell  were  able  to  correlate  transonic  drag  rise  data  for  peaky 
type  airfoils  on  older  Douglas  transport  aircraft  wings.  Their  correlation 
was  not  accurate  for  the  Boeing  747  or  the  F-111A/TACT  aircraft,  however. 

2.1  Empirical  Correlations  (C.R.  James) 

In  the  preliminary  design  phase,  gross  or  overall  correlations  are  used 
to  estimate  the  size  and  geometric  features  of  the  aircraft.  These 
correlations  are  traditionally  based  on  wind  tunnel  and  flight  test  data  from 
older,  successful  aircraft.  The  key  to  this  method  is  to  obtain  sufficient, 
relevant  data  from  aircraft  of  the  same  general  type  (bomber,  transport  or 
fighter,  etc.)  being  designed  in  order  to  form  trends  and  extrapolations  with 
confidence.  Data  collections  of  this  type  are  continuously  sought  and  main¬ 
tained  within  the  preliminary  design  groups  of  the  aircraft  industry.  The 
data  from  a  single  airframe  manufacturer,  for  a  particular  series  of  aircraft 
similar  to  the  new  design,  usually  produces  the  most  accurate  drag  prediction. 

Correlations  of  flight  test  data  based  on  aircraft  geometry  provide 
insight  as  to  the  relative  importance  of  factors  influencing  subsonic  minimum 
drag.  This  section  treats  two  types  of  aircraft:  bomber/transport  and 
fighter/attack.  Minimum  drag  is  related  to  wetted  area  through  the  use  of  an 
equivalent  skin  friction  coefficient,  C^g.  Fig.  7  compares  flight  derived 
values  of  C^  with  incompressible  flat  plate  turbulent  flow  skin  friction 
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coefficients  for  jet  aircraft  of  both  types.  This  comparison  shows  that,  in 
general,  about  two-thirds  of  subsonic  minimum  drag  may  be  attributed  to  skin 
friction.  The  balance  is  due  to  form  drag  and  interference.  Form  drag 
results  from  flow  over  curved  surfaces  and  flow  separation  induced  by  viscous 
effects.  Interference  drag  is  the  result  of  mutual  interaction  of  the  flow 
fields  developed  by  the  major  configuration  components.  The  equivalent  skin 
friction  coefficient,  C^e,  is  a  convenient  method  to  relate  total  drag 
(including  form  and  interference  drag)  to  wetted  area: 


q  ^wet 

where  D  =  Minimum  drag,  lb 

q  =  Dynamic  pressure,  Ib/sq  ft 

Swet  =  wettec*  area,  sq  ft 

The  quantity  D/q  is  defined  as  equivalent  parasite  area,  f. 

2.1.1  Bomber/Transport  Aircraft 

The  relationship  of  Cfg  and  f  to  wetted  area  is  illustrated  in  Fig.  8. 
Experimental  points  derived  from  flight  tests  of  several  transport  and  bombers 
are  indicated.  Wetted  areas  vary  from  approximately  5800  square  feet  for  the 
Breguet  941  to  33,000  square  feet  for  the  C-5A.  Equivalent  skin  friction 
coefficients  vary  from  0.0027  for  the  C-5A  to  0.0060  for  the  C-130.  The 
variation  in  Cfp  reflects  factors  which  influence  form  and  interference  drag 
as  well  as  Reynold's  number  effects  resulting  from  variations  in  aircraft 
size.  For  example,  the  three  turboprop  military  transports  (Br  941,  C-130, 
and  C-133)  exhibit  values  of  in  the  0.0050  to  0.0060  range.  These  air¬ 
craft  have  upswept  afterbodies  characteristic  of  rear  loading  arrangements  and 
consequently  incur  form  drag  penalties  which  increase  Cfg.  Interference 
effects  of  nacelles  and  propeller  slip  streams  also  contribute  to  the  higher 
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drag  levels.  Jet  transports  (737,  727,  and  DC-8),  with  a  higher  degree  of 
aerodynamic  cleanness,  are  in  the  0.0030  range.  The  C-5A,  because  of  its 
size,  operates  at  higher  Reynolds  numbers  and  has  the  lowest  of  the 
transports  evaluated.  The  low  of  the  C-5A  is  due  to  its  aerodynamic 
cleanness  and  very  low  aftbody  drag  increment,  despite  its  military  upsweep. 
2.1.1.  Fighter/Attack  Aircraft 

Correlation  of  flight  test  data  of  fighter/attack  aircraft  provides 
additional  insight  as  to  the  relative  importance  of  factors  influencing  sub¬ 
sonic  minimum  drag.  This  section  illustrates  effects  of  four  parameters: 
wetted  area,  relative  wing  size,  fuselage  shape,  and  aerodynamic  cleanness. 
Correlations  were  developed  from  flight  tests  of  eight  jet  fighters.  In  table 
2  the  aircraft  and  some  key  parameters  are  listed.  Wetted  areas  vary  from 
slightly  over  1100  square  feet  for  the  F-84G  to  approximately  2100  square  feet 
for  the  F-4E.  Equivalent  fuselage  fineness  ratios  vary  from  about  4.5  for  the 
F-86H  to  9.3  for  the  F-18.  These  aircraft  have  single  engines  with  the 
exception  of  the  F-4E/J,  F-14  and  F-18.  The  relationship  between  f,  C^g,  and 
wetted  area  is  illustrated  on  Fig.  9.  Experimental  points  derived  from  flight 
tests  of  the  aircraft  listed  on  Table  2  are  indicated  on  this  figure. 
Equivalent  skin  friction  coefficients  vary  from  0.0032  for  the  F-8D  to  0.0052 
for  the  F-4E/J.  The  variation  in  Cfe  results  primarily  from  factors  which 
influence  form  and  interference  drag.  Significant  factors  are:  (1)  relative 
wing  size,  (2)  fuselage  shape,  and  (3)  aerodynamic  cleanness. 

The  effect  of  relative  wing  size  is  illustrated  in  Fig.  10.  The  ratio  of 
the  equivalent  skin  friction  coefficient,  Cfe’  to  the  incompressible  mean 
flat  plate  value,  C^.  ,  is  shown  as  a  function  of  6,  the  ratio  of  the  total 

wetted  area  to  wing  wetted  area.  The  increase  in  this  ratio  with  increasing 
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TABLE  2 


AIRCRAFT  CHARACTERISTICS 
FIGHTER/ ATTACK  GEOMETRY 


MODEL 

DESIGNATION 

AREAS  (SQ. 

WING 

REFERENCE 

FT.) 

TOTAL 

WETTED 

EQUIVALENT 

FUSELAGE 

FINENESS 

RATIO 

A-7A 

375 

1691 

6.43 

F-4E/J 

530 

2092 

5.88 

F-8D 

375 

1821 

8.16 

F-84F 

325 

1257 

5.92 

F-84G 

260 

1104 

5.08 

F-86H 

313 

1186 

4.57 

F-100D 

400 

1571 

6.61 

F-105D 

385 

1907 

8.86 

F-18 

400 

2046 

9.28 

F-14 

565 

3097 

8.22 

14.5 


6  reflects  the  fact  that  wing  surface  tolerances  are  normally  held  much 
tighter  than  other  components,  and  protuberances,  roughness,  and  leakage 
penalties  are  smaller  for  the  wing  than  for  other  components. 

An  important  factor  affecting  fuselage  drag  is  the  degree  of  pressure  or 
form  drag  developed  due  to  flow  separation  induced  by  viscous  effects.  In 
general,  flow  separation  is  more  pronounced  with  bodies  of  low  fineness  ratio. 
This  effect  is  illustrated  in  Fig.  11  by  defining  an  equivalent  body  which 
includes  frontal  areas  of  the  fuselage,  canopy,  and  wing  to  establish  an 
equivalent  fuselage  fineness  ratio.  In  Fig.  12  this  parameter  is  illustrated 
and  defined.  A  form  factor,  Fg,  is  defined  from  the  empirical  expression: 

FB  =  1.02  1  +  1.5  +  _ 7 _ 

U/d)1*5  . 6 ( £./ d ) 3  (1-M3) 

Mach  number,  M,  accounts  for  compressibility  effects,  Fig.  11  also  includes  an 
aerodynamic  cleanness  factor,  w,  which  is  related  empirically  to  relative  wing 
size,  Fig.  13. 

The  aerodynamic  cleanness  factor  accounts  for  form  drag  and  interference 
effects  not  treated  by  fuselage  shape  parameters.  Examples  are  form  and 
interference  drag  for  wings  and  tails,  protuberance  and  leakage  drag,  manufac¬ 
turing  tolerance  effects,  and  surface  roughness.  The  magnitude  of  w  and  the 
variation  with  relative  wing  size  is  shown  in  Fig.  13.  Separate  trend  lines 
are  identified  as  a  function  of  the  number  of  inlet  and  exit  flows.  The  lower 
line  is  determined  by  the  single  inlet-single  nozzle  configurations,  the 
middle  line  by  the  data  from  the  F-105D  with  two  inlets  and  one  exit  nozzle, 
and  the  upper  line  by  the  data  from  the  F-4/J,  F- 14  and  F-18  with  two  inlets 
and  two  nozzles. 
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These  correlations  suggest  an  initial  approach  for  estimating  subsonic 
minimum  drag  for  similar  fighter/attack  configurations: 

1.  Calculate  the  wetted  area  and  a  characteristic  length  for  each 
component. 

2.  Calculate  an  area  weighted  mean  length,  1^,  for  the  complete  config¬ 
uration  or  estimate  1  from  experience. 

3.  Use  this  mean  length  to  calculate  Reynold's  number  at  the  flight 
condition  to  be  evaluated. 

4.  Determine  an  incompressible  skin  friction  coefficient,  C.^^,  using 
this  Reynold's  number.  The  Karman-Schoenherr  equation  for  turbulent  flow  or  a 
similar  expression  may  be  used. 

5.  Calculate  the  wetted  area  ratio,  6,  and  use  Fig.  13  to  establish  the 
value  of  a)  for  the  appropriate  inlet/nozzle  arrangement. 

6.  Calculate  the  equivalent  fuselage  fineness  ratio,  1/d,  and  determine 


the  form  factor,  Fg. 

7.  Determine  equivalent  parasite  area,  f,  from  the  expression: 
f  =  Cf  x  to  x  FB  x  Swgt 

ICO 

A  similar  approach  developed  for  bomber/transports  would  highlight  those 
factors  influencing  subsonic  minimum  drag  for  this  type  of  aircraft. 

2.2  Detailed  Drag  Estimates  -  Component  Buildup 


As  the  design  project  progresses,  the  geometric  definition  of  the  air¬ 
craft  is  formed  and  more  detailed  drag  prediction  methods  are  used  to  assess 
drag  targets  or  design  goals  for  each  aircraft  component.  These  goals  must  be 
consistent  with  the  overall  drag  target  previously  determined  since  the  total 
configuration  drag  is  the  sum  of  the  component  drags  including  interference 


effects. 
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Each  component  may  have  an  initial  run  of  laminar  flow,  followed  by  a 
transition  zone  and  fully  developed  turbulent  flow  thereafter.  The  extent  of 
each  flow  region  is  difficult  to  determine  with  precision.  Fully  turbulent 
flow  was  usually  assumed  from  the  leading  edge  aft  on  full-scale  flight 
vehicles  at  altitudes  below  about  twenty  kilometers.  The  smoothness  and 
rigidity  of  newer  composite  materials  may  provide  longer  runs  of  laminar  flow 
on  future  flight  vehicles.  Transition  is  usually  fixed  by  trip  strips  or  grit 
on  subscale  wind  tunnel  models  that  are  tested  at  lower  Reynolds  numbers.  The 
positioning  and  selection  of  boundary-layer  trip  strips  that  accurately 
simulate  full-scale  conditions  is  an  art.^’  ^ 

2.2.1  Friction  Drag 

There  exists  little  controversy  concerning  calculation  of  the  average 

laminar  skin  friction  on  one  side  of  a  doubly  infinite  flat  plate.  It  is 
g 

given  by  Blasius'  formula  from  his  exact  solution  to  the  laminar  boundary 
layer  equations  for  zero  pressure  gradient: 

Cf  =  1.32824/Re1/2 

where  ®  £  and  i  is  the  distance  from  the  stagnation  point  to  the 

v 

transition  zone. 

Transition  takes  place  when  the  length  Reynolds  number  nominally  exceeds 
3.5  x  10^  to  10^  for  a  flat  plate. ^  This  is  the  critical  Reynolds  number, 
often  assumed  to  be  one-half  million  for  many  flows,  although  pressure  gradi¬ 
ents  strongly  influence  the  location  of  transition  in  the  boundary  layer.  It 
is  the  objective  of  boundary-layer  stability  theory  to  predict  the  value  of 
the  critical  Reynolds  number  for  a  prescribed  main  flow.  Success  in  the 
calculation  of  the  critical  Reynolds  number  for  flows  of  general  aeronautical 
interest  has  eluded  workers  in  this  field  for  many  decades  despite  many 
dedicated  efforts. 
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Many  correlations  of  turbulent  boundary-layer  data  exist.  In  Fig.  14 
the  variation  of  with  Reynolds  number,  predicted  by  five  different  methods, 
is  displayed.  Significant  changes  in  can  occur  when  scaling  from  wind 
tunnel  (Rn  =  3  x  10®)  to  flight  (R  =  40  x  10®)  if  different  correlations  are 
used.  Two  well-known  and  widely-used  formulas  are: 

Prandtl-Schlicting: 

Cf  =  0.455  _  A_ 

(log  Re)2*58  Re 

where  A  depends  on  the  position  of  transition  and 
Karman-Schoenherr: 

Cf'1/2  =  4.13  log  (Re  Cf) 

The  Prandtl-Schlichting  formula,  being  explicit  in  Cf,  has  been  more 

108 

widely  used  in  the  past,  although  Paterson  has  shown  that  the  Karman- 

Schoenherr  formula  is  a  good  representation  of  existing  test  data  (Fig. 

15).  Both  formulas  have  been  used  as  the  basis  for  aircraft  drag  prediction 

methods.®8,  ^ A  new  method,  based  on  explicit  Prandtl-Schlichting  type 

153  154 

relations,  has  been  developed  by  White  and  Christoph  *  for  compressible 
turbulent  skin  friction.  Their  explicit  approximation 

Cf  =  0.42/ln2  (0.056  R  ) 

5  9 

is  accurate  to  ±  4%  in  the  Reynolds  number  range  10  -  10  .  It  is  the  formula 

122 

selected  by  Schemensky  for  his  complete  drag  prediction  method. 

The  reference  length  used  in  the  Reynolds  number  computation  is  the  body 

86 

length  for  near  bodies  of  revolution  and  MacWilkinson,  et.  al.  recommend  use 

of  the  mean  aerodynamic  chord  for  lifting  surfaces.  This  choice  of  reference 

length  will  be  less  accurate  (about  2%  low)  for  highly  swept  delta  wings.  A 

strip  method,  where  a  reference  length  is  calculated  for  each  streamwise  strip 

68 

of  the  wing,  should  be  used. 
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Compressibility  corrections  to  the  skin  friction  are  less  than  10%  at 
Mach  numbers  less  than  one  and  are  often  ignored  at  subsonic  speeds  unless 
continuity  with  supersonic  predictions  is  desired. 

2.2.2  Form  Drag 

Subsonic  minimum  drag  is  the  sum  of  the  friction,  form  and  interference 

drag  when  the  component  buildup  method  is  used.  The  form  drag,  or  pressure 

drag,  resulting  from  the  effects  of  nonzero  pressure  gradient  (component 

thickness)  is  usually  accounted  for  by  a  multiplicative  factor  applied  to  the 

122 

skin  friction  drag.  Following  Schemensky's  drag  accounting  system  for 
example, 

CD  Cf  *  Awet  *  FF  *  IF 

Components  Eq.  1 

+  CD  .  +  CD.  +  CD  . 

camber  base  mi  sc 

where  Awet  is  the  component  wetted  area 

FF  is  the  component  form  factor,  also  called  the  shape  factor,  SF 

IF  is  the  component  interference  factor 

The  interference,  camber,  base,  and  miscellaneous  drags  are  treated  in  later 

sections. 

The  prediction  of  the  form  drag,  or  FF,  is  dependent  on  empirical  corre¬ 
lations  as  given  by  Hoerner64,  Schemensky122,  DATCOM63,  O'Conner107, 
Snodgrass135,  MacWilkinson86,  and  others.  Example  formulas  for  the  compu¬ 
tation  of  FF  are: 

for  bodies:  FF  =  1  +  60/FR3  +  0.0025  *  FR 

for  nacelles:  FF  =  1  +  0.35/FR 
where  FR  =  Component  Length 
Width  *  Height 

and  for  airfoils: 
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FF  =  1  +  Kj  (t/c)  +  K2  (t/c)2  +  K3  (t/c)4 
where  t/c  is  the  streamwise  thickness  to  chord  ratio  and  the  1C  depend  on 
airfoil  series.  An  additional  term,  dependent  on  the  section  design  camber  is 
added  for  supercritical  airfoils. 

MacWilkinson,  et.al.  have  shown  that  fuselage  profile  drag  (the  sum  of 

skin  friction  drag  and  form  drag)  can  be  correlated  for  bodies  of  revolution 

169 

by  the  method  of  Young  ,  Fig.  16.  This  theory  agreed  with  experimental  data 

from  four  transport  fuselages  that  are  radical  departures  from  the  ideal  body 

of  revolution,  if  the  proper  fineness  ratio  is  used  as  shown  in  Fig.  17. 

93 

However,  excess  drag  due  to  upsweep  ,  which  may  amount  to  ten  percent  of  the 

total  cruise  drag  of  military  transport  aircraft,  is  not  apparent  from  these 

correlations.  A  new  prediction  method  for  reasonable  estimates  of  the  drag  of 

afterbodies  for  military  airlifters  has  recently  been  published  by  Kolesar  and 

May72  to  replace  the  method  of  Hahn,  et.al.50. 

Methods  of  estimating  airfoil  shape  factors  (or  form  factors)  were  also 

86 

compared  by  MacWilkinson,  et.al.  .  There  is  wide  scatter  in  both  the  theo¬ 
retical  and  experimental  values  of  the  shape  factor,  as  shown  in  Fig.  18, 
particularly  at  twelve  percent  thickness  to  chord  ratio,  a  nominal  value  for 
modern  transport  wings.  The  data  on  the  earlier  NACA  airfoils  were  not  suited 
to  accurate  assessment  of  airfoil  form  drag  because  of  the  testing  techniques 
used.  This  earlier  data,  however,  forms  the  basis  of  many  older  shape  factor 
correlations. 

Formulas  for  computing  the  form  factors  for  other  aircraft  components, 

64  1 99 

canopies,  stores,  struts,  etc.,  are  available  from  Hoerner  ,  Schemensky  , 
DATCOM53,  O'Conner*07,  and  Snodgrass*35,  along  with  examples  illustrating 
their  use  for  drag  estimation. 
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Formally,  the  form  drag  follows  the  same  dependence  on  Reynold's  number 

and  Mach  number  as  the  skin  friction  since  the  form  factors  are  constant  for  a 

given  geometry  in  most  drag  prediction  methods.  Hoerner^  has  developed 

compressible  form  factors  that  explicitly  account  for  Mach  number  effects  by 

transforming  the  appropriate  thickness  or  slenderness  ratio  using  the  Prandtl- 

2  1/2 

Glauert  transformation  factor,  (1  -  M  )  '  .  The  magnitude  of  these  correc¬ 
tions,  however,  precludes  verification  of  their  accuracy  due  to  lack  of 
precision  in  experimental  data. 

2.2.3  Interference  Drag 

Aerodynamic  interference  in  aircraft  is  the  change  in  flow  over  given 
components  of  the  configuration  due  to  the  presence  of  one  or  more  of  the 
other  components.  Interference  can  be  unfavorable  with  an  attendant  drag 
increase,  or  in  the  case  of  skillful  design,  favorable  where  the  sum  of 
component  drags  is  greater  than  the  total  drag  of  the  configuration.  The 
importance  of  interference  to  the  optimization  of  aircraft  to  meet  performance 
requirements  was  delineated  at  a  recent  AGARD  conference1^7.  Generally, 
important  interference  effects  involve  fluid  dynamic  phenomena  that  is  far  too 
complex  to  be  analyzed  by  existing  computational  methods. 

Interference  factors,  to  account  for  the  mutual  interference  between 
aircraft  components,  are  given  constants  based  on  experimental  evidence  for 
the  smaller  components  and  are  usually  presented  as  plotted  data  for  the  major 
components.  Representative  examples  of  this  can  be  found  in  Schemensky122  and 
DATCOM63. 

The  interference  factors,  unlike  the  form  factors,  are  correlated  as 
functions  of  Mach  and  Reynolds  numbers  for  the  major  aircraft  components. 

These  factors  are  usually  greater  than  unity  at  flight  Reynolds  number  since 
the  resultant  total  drag  is  usually  larger  than  the  sum  of  the  component  drags 
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when  tested  singly.  The  factors  also  increase  with  Mach  number  to  account  for 
compressibility  effects.  The  interference  factors  for  t.he  smaller  components 
are  constants! 

2.2.4  Camber  Drag 

This  increment  to  the  minimum  drag  accounts  for  camber  and  twist  in  the 
wing,  and,  in  some  drag  accounting  systems,  the  fact  that  all  aircraft  compo¬ 
nents  are  not  mounted  relative  to  each  other  to  attain  zero  lift  simultane¬ 
ously.  This  increment  is  also  related  to  the  drag  due  to  lift  methodology 


through  the  drag  accounting.  A  correlation  consistent  with  Equation  1  is 


°CAMBER 


_e_  K  (  C.  )2 
1-e  L 


where  e  =  l/(ir  *  AR  *  K)  span  efficiency  factor 


K  =  induced  drag  factor 

r  =  a  constant.  The  lift  coefficient  for  a  minimum  drag.  Often 
lL 

the  design  lift  coefficient. 

The  factors  e  and  K  are  determined  from  the  variation  of  drag  with  lift,  to  be 


considered  later. 


A  different  drag  predicting  method  has  been  shown  in  Fig.  2.  In  that 
system,  all  changes  in  profile  drag  with  lift  are  derived  from  experimental 
data  correlations  and  the  basic  profile  drag  is  not  incremented  for  camber 
directly.  The  method  of  Fig.  2  is  most  accurate  when  the  aircraft  designs  are 
closely  related  members  of  the  same  series,  for  example,  transports  from  a 
single  manufacturer. 

Thus,  drag  data  reduced  using  one  bookkeeping  system  may  not  include  all 
the  interference,  camber,  etc.,  effects  in  the  resulting  correlations  that  the 
drag  buildup  method  assumes  are  included,  unless  extreme  care  is  taken  to 
achieve  consistency. 
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2.2.5  Base  Drag 

This  term  is  a  weak  function  of  Mach  number  at  low  speed  and  is  given  by 

Lj  64 
Hoerner  as: 

VsE  =  (0'1  +  °'1222H8)  SBASE/SREF 

It  obviously  gains  significance  with  increasing  Mach  number,  however,  it  is 
independent  of  lift  and  is  often  not  correlated  separately  from  the  body  drag. 

2.2.6  Miscellaneous  Drag 

This  term  accounts  for  the  differences  between  the  corrected  aerodynamic 

reference  model  drag  and  the  full-scale  airplane  drag  due  to  surface  irregu- 

larities,  such  as  gaps,  mismatches,  fasteners  and  small  protuberances  ,  and 

leakage  due  to  pressurization.  It  is  estimated,  from  experience  and  aircraft 

type,  as  some  percentage  of  the  total  friction,  form  and  interference  drags  in 

the  preliminary  design  stages.  In  the  later  design  stages,  when  configuration 

details  are  known,  each  of  these  small  drag  terms  may  be  accounted  for  sepa- 

170 

rately.  An  AGARD  conference  ,  devoted  to  this  subject,  was  held  in  1981. 

The  transfer  of  wind  tunnel  data  to  full-scale  flight  conditions  is  discussed 
under  wind  tunnel/flight  correlation  of  lift,  drag,  and  pitching  moment. 

It  would  be  highly  desirable  to  assess  the  creditability  of  each  term  in 
the  drag  buildup  equations.  Unfortunately,  for  complete  aircraft  configura¬ 
tions,  this  would  require  the  ability  to  determine  accurately  the  variation, 
with  Reynolds  and  Mach  numbers,  of  the  drag  forces  due  to  viscous  and  pressure 
effects  of  each  aircraft  component  -  an  intractable  task. 

2.2.7  Drag  Due  to  Lift 

Various  methods  have  been  used  to  correlate  subsonic  drag  due  to  lift,  or 
induced  drag,  depending  upon:  the  design  project  phase,  the  configuration, 
the  anticipated  range  of  angle  of  attack,  and  the  sophistication  of  the 
prediction  method. 
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The  complexity  of  drag  due  tc  lift  prediction  methodology  is  illustrated 

lOO 

in  Fig.  19,  from  Schemensky  L,  based  on  the  extensive  systematic  correlations 
131 

of  Simon  .  The  drag  is  predicted  in  each  of  the  regions,  one  through  seven, 

by  modified  or  distinct  equations  and,  perhaps,  additional  correlations 

(terms)  to  account  for  increasing  complexity  in  the  flow  field.  The  following 

section  describes  drag  calculation  in  regions  one,  five,  and  six  below  the 

critical  Mach  number,  M^.  The  detailed  drag  methodology  is  contained  in 
122 

Schemensky 's  report. 

Classical  aerodynamic  theory  predicts  that  the  drag  due  to  lift  is  a 

function  of  the  vortex  distribution  shed  behind  the  wing  (downwash).  It  is 

parabolic  with  lift  and  is  given  by 

C  =  C2  h  AR 
ui 

for  an  elliptic  lift  distribution  spanwise.  Wing  taper  and  sweep  produce 

rcrielliptic  loadings,  thus 

Cn  =  C2  /ti  AR  e' 
ui  1 

where  e'  <  1.0  is  the  wing  alone  efficiency  factor.  An  empirical  modification 
to  e‘  to  account  for  the  presence  of  the  fuselage  in  the  spanwise  lift 


distribution  results  in 


C2  /-rr  AR  e  =  N  C2 


where 


e  =  e1  [1.  -  (d/b)  ] 


and  d/b  is  the  body  diameter  to  wingspan  ratio.  The  effects  of  camber  are  not 

included  in  these  equations.  Camber  causes  a  shift  in  the  polar  so  that 

minimum  drag  occurs  at  a  finite  lift  coefficient,  or  aC^  depending  on  the 

131  0 

author  of  the  method.  Simon,  et.al.  illustrated  this  camber  effect  and  the 
nonparabolic  nature  of  the  "drag  polar"  at  higher  lift  coefficients  where  flow 


separation  effects  become  important,  as  shown  in  Fig.  2C. 
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Simon  found  a  further  empirical  modification  to  this  equation  was 

necessary  to  account  for  local  supervelocities  at  the  wing  leading  edge,  while 
122 

Schemensky  ,  using  the  camber  drag  term,  found  that 

K  =  L-R  +  R 

irAR  e 

Here  R  is  the  leading-edge  suction  parameter  and  properly  bounds  K  between 

1/irAR  eo  for  full  leading-edge  suction  (R  =  l.)»  and  the  upper  bound  of  drag, 

1/C^,  for  zero  leading-edge  suction,  or  tan  a.  The  suction 

parameter,  R,  is  significantly  affected  by  airfoil  camber,  conical  camber, 

leading-edge  radius,  Reynolds  number,  and  sweep.  The  proper  determination  of 

122 

R  is  an  involved  procedure. 

It  is  within  region  one,  bounded  by  the  onset  of  leading-edge  separation 

with  reattachment,  i.e.,  the  for  initial  polar  break  boundary,  and  the 

critical  Mach  number  boundary  that  surface  paneling  methods  are  accurate  and 

can  be  applied.  These  methods  are  all  based  on  the  Prandtl-Glauert  equation, 

a  linearized  form  of  the  complete  fluid  dynamics  equations  of  motion,  and  can 

predict  the  lift  coefficient  for  minimum  drag  and  the  variation  of  drag  with 

lift  for  complex,  complete  aircraft  configurations  quite  accurately. 

Four  panel  methods  were  systematically  evaluated  in  1976  at  an  AGARD 
140 

meeting  ;  however,  many  new  methods  are  currently  in  use;  for  example: 
DOUGLAS-HESS58’  59,  MCAERO16,  PAN  AIR26,  39’  130,  146,  QUADPAN171,  USSAERO41, 
and  VSAERO38’  88.  Each  has  its  particular  advantages  and  limitations^43, 
in  addition  to  those  due  to  the  Prandtl-Glauert  equation,  because  of  the 
particular  numerical  method  used.  An  evaluation  of  these  six  newer  methods 
against  the  datum  cases  of  Reference  140  would  be  enlightening.  Drag  predic¬ 
tion,  not  specifically  considered  during  that  meeting,  should  be  added  as  an 
evaluation  criteria. 
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Basically,  all  newer  panel  methods  use  Green's  theorem  and  assume  an 

algebraic  form  for  the  source  and  doublet  singularity  strengths  on  each  of  the 

panels  that  are  used  to  represent  the  surface  of  the  aircraft.  Higher-order 
16  26  59 

methods  5  ’  use  polynomials  to  represent  the  singularities  while  lower- 

order  methods®®’  171  assume  that  the  singularities  are  constant  (or  vary 
linearly  in  one  direction)  over  each  panel.  Even  though  the  accuracy  of  each 
type  of  representation  increases  as  the  paneling  density  is  increased,  so  does 
the  computer  time  and  cost,  thus  forcing  the  usual  trade-off  between  accuracy, 
schedule  and  cost. 

Panel  methods  have  primarily  been  used  as  flow  diagnostic  tools®* ■*  by 
performing  systematic  comparative  studies  of  alternative  geometries  and  noting 
regions  of  high  velocity  or  velocity  gradient  for  further  detailed  analysis. 
The  proper  application  of  any  panel  method  remains  an  art  due  to  the 

nc 

myriad  of  paneling  layouts  that  can  be  used  to  represent  the  configurations  ’ 

,  the  impact  of  various  types  of  boundary  conditions  ’  and  the  many 

59  143  146  171 

numerical  ways  of  enforcing  the  Kutta  condition  ’  *  ’  The  calcu¬ 

lation  of  the  induced  drag  from  the  velocity  potential  also  presents  problems. 
The  near-field  summation  of  pressure  multiplied  by  the  area  projected  in  the 
drag  direction  results  in  differences  of  large  numbers  and  a  loss  of  signifi- 
cant  figures.  Control  volume  and  the  associated  far-field  or  Trefftz  plane 
wake  methods®®’  7^  for  drag  calculations  also  present  conceptual  and  numerical 
problems  when  applied  to  complete  wing,  fuselage  and  tailplane  configurations. 
The  optimum  panel  method  primarily  depends  upon  the  skill  of  the  user,  the 
computational  resources  available,  the  geometric  configuration  to  be  analyzed, 
and  the  type  and  accuracy  of  the  results  desired. 


26 


Although,  drag  prediction  results  are  not  often  published  in  the  open 
literature,  two  examples  have  been  found  to  demonstrate  the  induced  drag 
prediction  capabilities  of  panel  methods. 

98 

The  first,  from  Miller  and  Youngblood  ,  is  shown  in  Fig.  21.  The 
intricate  paneling  required  for  the  advanced  fighter  configuration  is  shown, 

ion 

along  with  the  drag  polars  predicted  using  PAN  AIR  °  and  a  lower-order 
165 

method,  USSAERO-B  .  The  minimum  drag  was  predicted  using  the  empirical 
methods  of  the  previous  sections.  The  comparison  with  experimental  data  is 
good,  with  the  results  favoring  use  of  the  higher-order  PAN  AIR  method; 
however,  the  computer  resources  necessary  to  achieve  the  increased  accuracy 
were  considerably  higher.  Modeling  studies,  particularly  with  respect  to  the 

no 

placement  of  the  canard  wake,  were  also  necessary  to  achieve  these  results  . 

31 

The  second  example  is  from  Chen  and  Tinoco  .  The  versatility  of  the  PAN 

AIR  method  in  predicting  the  effects  of  engine  power  on  the  lift  and  drag  of  a 

transport  aircraft  wing,  body,  strut,  and  nacelle  wind  tunnel  model  is  shown 

in  Fig.  22.  The  computed  lift  and  drag  increments  on  various  components  of 

the  aircraft  due  to  increasing  the  fan  nozzle  pressure  ratio  from  "ram"  to 

"cruise"  shows  that  the  wing,  the  fuselage,  and  the  strut  all  had  favorable 

contributions,  while  the  nacelle  was  the  only  component  that  contributed  to  a 

drag  penalty.  When  the  lost  lift  was  restored  and  the  associated  induced  drag 

was  included,  the  computed  blowing  drag  was  1.3  drag  counts.  About  two  drag 

counts  were  measured  in  the  wind  tunnel. 

Additional  applications  of  panel  methods  include:  fighter  aircraft  with 

28 

externally  carried  stores  (mutual  interference  effects),  transport  aircraft 

102 

with  flaps  deployed  ,  and  separated  flow  modeling  using  free  vortex 
88 

sheets  .  Panel  methods  have  also  been  used  with  boundary-layer  calculation 
methods31,  70  to  simulate  viscous  displacement  effects  and  the  parabolized 
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Navier-Stokes  equations'^.  These  hybrid  methods,  designed  to  conserve 
computer  resources  by  using  only  the  equations  necessary  to  analyze  the  flow 
in  each  appropriate  zone,  have  had  limited  success  in  predicting  the  fluid 
velocities  and  pressure^'  and  lift  loss  due  to  viscous  effects'^.  Drag 
prediction  is  presently  well  beyond  their  grasp. 

Drag  due  to  lift  prediction  in  regions  five  and  six  of  Fig.  19  is  an 
empirical  art  and  it  is  likely  to  be  many  years  before  satisfactory  physical 
and  mathematical  models  are  achieved  for  the  complex  separated  flows  that 
occur  in  these  regions  For  example,  the  lift  coefficient  for  polar  break, 
CLr>  ,  and  the  initial-stall  lift  coefficient,  CLd  ,  coincide  for  thin  wings, 

*  Ij 
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low  Reynolds  number,  or  highly  swept  wings,  and  region  five  does  not  exist 

The  lower  bound  for  region  five  is  the  lift  coefficient  at  which  leading- 
edge  separation  and  subsequent  reattachment  occur  -  the  bubble-type  separa¬ 
tions.  The  upper  bound  is  reached  when  stall  occurs  and  the  entire  upper 
surface  flow  is  separated.  Blunt,  thick  airfoils  generally  exhibit  trai ling- 
edge  separation,  while  very  thin  airfoils  exhibit  leading-edge  separation. 
Airfoils  of  moderate  thickness  are  likely  to  separate  and  reattach  at  the 
leading-edge,  followed  by  trailing-edge  separation  and  stall  at  higher  lift 
coefficients.  The  leading-edge  bubble  separation  produces  an  increase  in  drag 
due  to  lift  because  of  decreased  leading-edge  suction.  Above  drag  break,  the 

flow  separates  completely  from  the  wing  and  the  drag  rapidly  increases. 

122  131 

The  drag  correlations  developed  by  Schemensky  11 ,  Simon,  et.al.  ,  and 
Axelson^  can  be  used  for  rough  estimates  during  the  preliminary  design  phase, 
but  lack  the  logical  and  consistent  framework  necessary  for  confident,  de¬ 
tailed  design  trade-off  studies.  As  shown  in  Fig.  20,  the  drag  due  to  lift  in 
region  five  varies  with  lift  at  a  greater  rate  than  in  region  one.  The 
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additional  drag  in  region  six  is  a  separation  drag  increment  that  has  been 

131 

correlated  with  the  lift  coefficient  at  drag  break 

Numerous  symposia  have  been  devoted  to  the  prediction  of  viscous-inviscid 

45 

interactions,  separation  and  stall.  Gilmer  and  Bristow  accurately  predicted 

the  lift  and  drag  of  three  airfoils  through  stall  by  combining  a  two- 

dimensional  panel  method,  an  integral  boundary-layer  method  and  an  empirical 

relationship  for  the  pressure  in  the  separated  zone  and  wake  recovery  region. 

89 

Maskew,  Rao,  and  Dvorak  ,  also  using  a  panel  method  with  a  simple  wake  model 
with  prescribed  separation  line  and  wake  geometry,  were  able  to  predict  the 
pressure  distribution  over  an  aspect  ratio  six,  10°  swept  wing  at  21c  angle  of 
attack. 

The  subject  of  leading-edge  separation  from  slender  wings  has  recently 

rn 

been  summarized  by  Hitzel  and  Schmidt  .  Potential  flow  methods  that  model 

the  leading-edge  vortex  and  recently-developed  Euler  methods  can  predict  the 

flowfield  and  lift,  but  drag  comparisons  are  not  included.  The  Pohlhamus 

24 

suction  analogy  in  conjunction  with  linear  lifting  surface  theory  has 
produced  excellent  comparisons  with  experimental  drag  data  for  thin  slender 
wings.  However,  the  method  fails  for  higher  aspect  ratio  wings  typical  of 
subsonic  transport  aircraft. 

The  use  of  strakes  and  other  leading-edge  devices  to  reduce  the  drag  at 
high  lift  and  transonic  speeds  typical  of  maneuvering  fighter  aircraft  was  the 
subject  of  an  extensive  AGARD  conference61.  Design  guidelines  and  empirical 
rules  for  strakes  were  determined  from  extensive  wind  tunnel  data  bases.  A 
fundamental  understanding  of  the  effects  of  strakes  on  the  wing  flowfield, 
lift  ano  drag  was  not  evident  in  the  proceedings. 

High  angle  of  attack  aerodynamics61  is  certainly  one  area  of  aerodynamics 
that  is  not  a  mature  science.  Considerable  theoretical  and  experimental 
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effort  will  be  required  to  develop  accurate  force  and  moment  prediction 
methods  for  general  configurations. 

3  Transonic  Drag 

The  format  of  this  section  parallels  that  of  Section  2.  A  rapid, 
strictly  empirical  method  for  determining  the  overall  drag  rise  for  prelimi¬ 
nary  design  is  presented,  followed  by  a  drag  component  buildup  method  to 
assess  drag  targets  for  each  aircraft  component.  The  section  concludes  with  a 
summary  of  the  status  of  numerical  aerodynamics  methods  for  detailed  drag 
cleanup  during  the  later  design  phases. 

3.1  The  Drag  Rise  (C.R.  James) 

A  correlation  of  the  point  drag  difference  between  Mach  numbers  1.2  and 
0.8  for  many  fighter  aircraft  is  shown  in  Fig.  23.  The  shape  of  the  drag 
curve  between  these  Mach  numbers,  where  the  drag  may  attain  a  maximum,  cannot 
be  determined  by  this  method. 

The  increase  in  minimum  drag  coefficient  at  transonic  speeds  is  deter¬ 
mined  primarily  by  configuration  slenderness.  Body  fineness  ratio  (length/- 
diameter),  lifting  surface  sweeps,  and  thickness  ratio  (thickness/chord 
length)  are  primary  parameters  for  expressing  slenderness  of  major  configura¬ 
tion  elements.  The  area  rule  concept^’  ^  provides  a  convenient  approach  to 
combining  body  and  wing  cross-sectional  areas  to  define  an  equivalent  body  of 
revolution.  Fineness  ratio  of  this  equivalent  body  is  used  as  an  index  of 
configuration  slenderness  to  derive  a  correlation  of  transonic  drag  rise. 

Historic  data  from  twenty-eight  aircraft,  where  the  cross-sectional  area 
distribution  normal  to  the  body  centerline  (Mach  1.0  cut)  is  used  to  define  an 
equivalent  body  fineness  ratio  U/d),  is  compared  in  Fig.  23.  In  all  air- 
breathing  configurations,  the  inlet  capture  area  is  subtracted  from  the  body 
total  area  aft  of  the  inlet  station  to  define  the  equivalent  body  shape. 
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Minimum  drag  coefficient,  C^,  is  based  on  the  maximum  cross-sectional  area  of 
the  equivalent  body  and  the  difference  between  drag  values  at  M  =  1.20  and  M  = 
0.80  is  defined  here  as  the  drag  rise. 

Three  generations  of  aircraft  are  shown:  (1)  early  experimental  air¬ 
craft,  (2)  the  Century  series  aircraft,  and  (3)  subsequent  developments.  Most 
of  the  data  are  taken  from  flight  test  results,  although  data  for  some  air¬ 
craft  (the  NASA  LFAX  configurations  and  the  XFY-1,  for  example)  are  taken  from 
wind  tunnel  tests.  The  theoretical  slender  body  values  determined  for  Sears- 

1  nc 

Haack  minimum  drag  body  shapes,  pointed  on  both  ends,  is  indicated  as  a 
lower  limit.  An  overall  progression  toward  this  lower  limit  is  discernable. 
The  first  generation  supersonic  aircraft  with  fineness  ratios  of  the  order  of 
7.0  had  drag  rise  values  about  60%  higher  than  the  Sears-Haack  limit,  the 
second  generation  about  50%,  and  the  third,  about  25%.  The  B-58  with  an 
external  store  pod  is  the  configuration  nearest  the  lower  limit. 

Unfortunately,  some  of  the  currently  operational  aircraft  depart  from  the 
trends  established  by  the  earlier  generations.  The  F-15,  for  example,  with  a 
fineness  ratio  of  8.7  has  a  drag  rise  200%  higher  than  the  Sears-Haack  lower 
limit  and  the  F-18  with  a  fineness  ratio  of  9.5  is  225%  above  this  limit. 

3.2  Detailed  Drag  Estimates  -  Component  Buildup 
3.2.1  Zero-Lift  Drag 

A  schematic  diagram  of  the  component  buildup  method  is  shown  as  Fig. 

24.  The  skin  friction  drag  may  be  computed  by  the  methods  of  Section  2.2  at 
the  appropriate  Reynolds  and  Mach  numbers  for  the  flow.  The  friction  drag  is 
nearly  invariant  over  the  small  range  of  Mach  number  that  is  considered 
transonic,  unless  the  localized  details  of  shock  boundary-layer  interactions 
are  considered. 
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While  the  form  and  interference  drags  are  also  computed  by  the  methods  of 
Section  2.2,  both  are  assumed  to  terminate  abruptly  at  sonic  speed.  Recall 
that  some  form  factors  are  constants  over  the  entire  subsonic  range  while 
others  are  functions  of  both  Reynolds  and  Mach  numbers.  Similarly,  the  wave 
drag,  to  be  treated  in  Section  4.2,  is  assumed  to  begin  abruptly  as  Mach  one 
is  exceeded.  Since  neither  assumption  is  physically  realistic,  the  zero-lift 
drag  curve  is  forced  to  be  continuous  through  the  first  derivative  by  fitting 
a  polynomial  through  the  zero-lift  drag  value  at  critical,  or  drag  divergence, 
Mach  number.  Dreg  creep  is  accounted  for  by  the  variation  of  form  and 
interference  factors  with  subsonic  Mach  number  given  in  Section  2.2.  Thus, 
the  drag  rise  and  the  interference  plus  form  drag  are  replaced  by  wave  drag 
beyond  Mach  one. 

The  decrease  in  critical  Mach  number  with  increasing  lift  for  most 
airfoils  causes  the  subsonic  drag  polar  to  begin  its  rise  earlier  at  higher 
lift  coefficients.  The  drag  rise  is  usually  separated  into  a  minimum  drag 
contribution  and  a  contribution  due  to  lift  to  account  for  this  change.  The 
drag  rise  due  to  lifting  surfaces  begins  at  MCr  while  the  drag  rise  due  to  all 

ether  components  begins  at  Mr  as  correlated  by  Schemenskv1'1  . 

Ro 

Certainly,  many  unresolved  problems  and  contradictions  are  evident  in 
this  method.  It  is  merely  an  acceptable  artifice  to  provide  continuous  drag 
curves  through  the  transonic  regime.  For  example,  the  component  critical  Mach 
numbers  depend  on  airfoil  section,  thickness  ratio,  or  body  slenderness 

•?  o  o 

ratio1  ,  It  may  be  different  for  each  aircraft  component,  causing  each  to 
enter  the  transonic  regime  at.  a  different  free  steam  Mach  number,  or  interfer¬ 
ence  between  closely  spaced  components,  such  as  aft  pylon-mounted  engines,  may 
dominate  the  entire  flow  field  and  drag'"''7.  The  drag  curves  generated  by  this 
curve  fitting  also  do  not  account,  for  the  pressure  drag  slope  relationships 
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derived  from  transonic  similarity  theory  ,  and  so  on.  Transonic  drag 
prediction  is  not  yet  a  reliable,  consistent  art! 

3.2.2  Drag  Due  to  Lift 

In  the  transonic  region  three  of  Fig.  19,  bounded  by  the  critical  Mach 
number  and  the  limit  Mach  number  ,  where  0.95  <  <  1.00,  the  induced 

drag  is  computed  in  the  same  manner  as  described  for  region  one,  Section  2.2. 
The  drag  rise  for  lifting  surfaces  from  the  previous  section  shifts  the  basic 
polar,  depending  upon  MCr  .  In  the  transonic  region  four  above  M^,  the  drag 
polar  is  calculated  by  interpolation  of  the  polar  shape  factors,  K,  in  regions 
two  and  three. 

A  newer  method  for  transonic  drag  due  to  lift  prediction  to  extremely 

4 

high  angles  of  attack  has  been  programmed  by  Axelson  ,  This  method  is  based 
on  new  forms  of  compressible  wing  theory  covering  potential  and  nonpotential 
flows  and  requires,  in  addition  to  the  usual  geometric  inputs,  the  limiting 
forward  and  aft  chordwise  locations  of  the  shock  wave.  The  method  predicts 
the  drag  due  to  lift  reasonably  well  for  quite  general,  assumed  limit  shock 
positions,  as  shown  in  Fig.  25.  The  drag  at  zero  lift  is  not  calculated  by 
this  program  and  must  be  supplied  by  the  user. 

3.3  Numerical  Transonic  Aerodynamics 

The  progress  in  computational  transonic  aerodynamics  has  been  reviewed  in 
several  recent  publications106’  167  ’  168.  The  book,  edited  by  Nixon106, 

describes  steady  transonic  flow  research  from  both  an  industry  viewpoint  and  a 
research  viewpoint  and  includes  recent  advances  in  experimental  techniques  and 
prediction  methods.  Yoshihara's  AGARD  reports167’  168  review  transonic  aero¬ 
dynamics  for  the  applied  aerodynamicist  engaged  in  the  design  of  combat  and 
airlift  aircraft  with  proper  emphasis  on  aircraft  component  interference167. 
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The  numerous  authors  in  the  field  of  numerical  transonic  aerodynamics 
have  reached  a  consensus  —  transonic  drag  predictions  are  currently  unreli¬ 
able  by  any  method.  Neither  the  polar  shape  nor  drag  levels  are  acceptably 
predicted.  The  uncertainty  in  the  magnitude  and  often,  even  the  sign,  of 
incremental  drag  due  to  small  geometric  changes  precludes  the  use  of  drag  as 
an  object  function  in  minimization  procedures.  The  surrogate  transonic  design 
criteria  is  the  pressure  distribution  which  places  a  heavy  burden  on  the 
aerodynamicist.  The  target  pressure  distribution  must  be  chosen  such  that 
boundary-layer  separation  is  avoided,  wave  and  vortex  drag  are  minimized,  and 
the  resulting  geometry  is  acceptable. 

Transonic  flow  is  highly  three-dimensional  and  inherently  nonlinear.  A 
heirarchy  of  approximate  equations,  in  addition  to  the  full  Navier-Stokes 

equations,  have  been  used  to  analyze  these  flows  with  varying  degrees  of 
168 

success  .  Unfortunately,  the  geometric  fidelity  of  the  numerical  represen¬ 
tation,  necessary  to  distinguish  drag  differences  due  to  design  changes, 
diminishes  as  the  higher  approximations  are  used.  This  is  a  result  of 
increased  grid  generation  and  computational  storage  requirements.  Grid 

resolution  has  been  a  primary  cause  for  poor  drag  prediction  using  even  the 
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full  Navier-Stokes  equations  with  various  turbulence  models  . 

The  theoretical  and  numerical  bases  for  drag  calculations  based  on  the 

approximate  equations  (transonic  small  disturbance,  full  potential  and  Euler) 

have  been  questioned1  *  1  *  c  and  remain  controversial.  The  addition  of  a 

coupled  boundary-layer  solver  to  the  inviscid  transonic  flow  calculation 

1 68 

further  complicates  the  problem  of  accurate  drag  determination  . 

Transonic  flow  calculations  are  presently  limited  in  geometric  capability 
to  various  collections  of  aircraft  components  and  lack  the  full  configuration 
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geometric  capability  of  panel  methods.  The  following  examples  indicate  the 

current  limitations  and  progress  that  is  being  made  in  this  area. 

84 

Longo,  Schmidt,  and  Jameson  illustrated  the  improvement  in  transonic 
airfoil  drag  prediction  produced  by  modifying  the  original  Bauer-Garabedian- 
Korn-Jameson  (BGKJ)  nonconservative  full  potential  method  with  weak  viscous 
interaction  to  account  for  strong  viscous  interactions  (VGK)  as  shown  in  Fig. 
26.  The  Dofoil  method  consists  of  a  conservative,  full -potential  flow  solver 
coupled  to  a  set  of  integral  boundary-layer  methods  with  special  models  for 
separation  bubbles  at  the  trailing  edge.  Further  improvements  (unpublished) 
have  been  obtained  by  substituting  a  finite  volume  Euler  flow  solver  for  the 
inviscid  method. 

57 

Henne,  Dahlin  and  Peavey  ,  also  using  the  Bauer,  Garabedian  and  Korn 

method,  obtained  an  excellent  correlation  of  drag  divergence  Mach  number. 

They  found  significant  scatter  in  the  drag  level  and  compressibility  drag 

increment  results  for  eight  airfoils. 

60 

Hicks  extensive  review  of  transonic  wing  design  using  both  conservative 

and  nonconservative  potential  flow  codes  did  not  use  drag  as  an  evaluation 

criteria  since  it  was  considered  to  be  the  least  accurate  quantity  calculated 

by  the  aerodynamic  codes.  Poor  experiment  to  theory  correlations  of  the  wing 

pressure  distributions  were  found  for  configurations  with  low  fineness  ratio 

57 

bodies.  Henne,  Dahlin,  and  Peavey  again  found  drag  divergence  Mach  number 
could  be  correlated  while  drag  level  and  compressibility  drag  increment  could 
not,  for  three  high  aspect  ratio  transport-type  wings,  as  shown  in  Fi  9.  27. 

i^o 

Waggoner  ,  using  a  small  disturbance  transonic  analysis  code  coupled  with  a 
two-dimensional  boundary-layer  code  for  the  analysis  of  a  supercritical 
transport  wing-body  configuration  found  the  method  sensitive  to  design 
changes,  but  an  increment  of  41  drag  counts  was  added  to  the  predicted  drag 
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level  to  obtain  agreement  with  test  data.  The  computational  grid  resolution 
was  coarse,  particularly  or,  the  body  nose  and  boattail  region  which  may 
account  for  this  drag  increment. 

The  preceding  results  show  that  transonic  aerodynamics  has  enjoyed 
limited  success  when  applied  to  high  aspect  ratio  transport  configurations. 
Hov/ever,  when  applied  to  a  high-performance  fighter  configuration  with  a  wide 
complex  fuselage,  or  canards,  none  of  the  ful 1 -potential  methods  is  capable  of 
straight-forward  prediction  of  the  wing  pressures,  even  at  the  outboard  span 
stations. 

With  the  possible  exception  of  the  development  of  a  new  transonic  method 
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based  on  paneling  technology  ,  geometric  generality  has  been  limited  to  the 

127 

small  disturbance  equations.  Shanker  and  Goebel  have  developed  a  numerical 
transformation  capable  of  modeling  the  closely-coupled  canard-wing  inter¬ 
actions  typical  of  proposed  highly  maneuverable  fighter  configurations.  Drag 

90 

minimization  for  these  complex  configurations  has  been  addressed  by  Kason 
using  a  modified  vortex  lattice  method. 

Transonic  aerodynamic  theory  and  the  associated  numerical  methods  are 
currently  employed  for  design  guidance.  It  will  be  some  time  until  transonic 
drag  predictions  are  accurate  and  reliable  enough  for  detailed  design  data. 

4  Supersonic  Drag 

Somewhat  contrary  to  the  state  of  transonic  drag  prediction,  supersonic 
drag  is  firmly  based  on  theoretical  and  numerical  methods.  In  the  early 
1970s,  the  supersonic  transport  and  the  B-70  programs  provided  the  impetus  for 
the  extensive  development  and  numerical  implementation  of  linearized  super¬ 
sonic  theory  that  had  begun  in  the  1960s.  The  dedicated,  pioneering  efforts 
of  Bonner  5  *  *  ,  Carlson^  ,  Carmichae r " ,  Harris'  5  ,  Middleton^  , 


36 


1  f>5 

Woodward  '  ,  and  others  have  produced  automated  numerical  design  procedures 
comprised  of  many  individual  computer  codes  under  the  control  of  an  executive 
coupling  routine.  These  numerical  procedures  are  continually  evolving  to 
reduce  the  shortcomings  uncovered  during  the  succeeding  years  of  intensive 
application  to  fighter,  bomber  and  transport  configurations.  The  drag  pre¬ 
dictions  produced  by  the  current  versions  of  these  programs  are  excellent  for 
many  slender  supersonic  configurations  at  low  lift  coefficients  and  are  widely 
used  throughout  the  aircraft  industry  for  preliminary  design. 

As  shown  in  Fig.  5a,  the  supersonic  cruise  drag  components  are  dramat¬ 
ically  different  from  the  subsonic  cruise  drag  components.  The  large,  obvious 
difference  is  the  wave  drag,  a  component  that  does  not  have  a  counterpart  in 
subsonic  flow.  The  lift-dependent  and  skin  friction  drag  components  are 
analogous  to  their  subsonic  flow  counterparts  and  are  calculated  by  the 
following  modifications  to  the  previous  subsonic  methods.  The  supersonic  drag 
buildup,  using  linearized  theory,  is  illustrated  in  Fig.  28. 

4.1  Friction  Drag 

The  effects  of  compressibility  and  heat  transfer  must  be  accounted  for  at 
supersonic  speeds.  The  mean  adiabatic  skin  friction  coefficient  decreases  by 
about  30%,  regardless  of  Reynolds  number,  in  changing  from  incompressible 
speeds  to  Mach  two,  as  shown  in  Fig.  29.  The  relative  predictions  from  the 
four  theories  vary  with  both  Reynolds  and  Mach  numbers,  illustrating  several 
drag  prediction  problems.  Certainly,  the  basic  full-scale  skin  friction  drag 
prediction  during  initial  design  depends  on  the  theory  selected.  The  corre¬ 
lation  from  wind  tunnel  model  scale  to  full-scale  flight,  as  illustrated  by 

O] 

the  Concorde  data,  depends  on  the  theory  selected.  Finally,  the  reduction 
of  flight  test  results,  in  order  to  define  actual  drag  as  a  function  of  Mach 
number,  altitude,  ambient  temperature  and  lift  coefficient,  is  also  theory 


37 


dependent.  The  magnitude  of  the  drag  change  produced  by  selection  of  alter- 

p  1 

nate  prediction  methods  is,  however,  small.  Leyman  and  Markham  have  shown 
that  compared  with  the  later  method  of  Winter  and  Gaudet,  Michel's  method 
predicted  lower  aircraft  drag  directly  by  2-1/2?.  C,  (0.8%  total  drag)  and 
lower  aircraft  drag  from  extrapolation  of  the  l/45th  scale  model  data  by  5%  Cp 
(1.5%  total  drag)  for  the  Concorde.  The  Winter  and  Gaudet  method  is  similar 
to  that  of  Sommer  and  Short10  ,  and  each  reproduced  available,  acceptable, 
experimental  results  more  accurately  than  Michel's  method.  However,  the 
scatter  in  all  data  precluded  changing  skin  friction  prediction  methods  during 
the  Concorde  design  process.  The  same  theory  should  be  used  for  all  estimates 
during  the  design  process  to  achieve  consistent  results. 

The  Sommer  and  Short  method  is  used  to  calculate  the  compressible  skin 
friction  coefficient,  Cp,  from  a  reference  skin  friction  coefficient,  Cp,  for 
a  selected  free  stream  Mach  number,  M^,  and  Reynolds  number,  R^,  and  adiabatic 
wall  temperature,  Tw,  using  the  incompressible  Karman-Schoeriherr  formula 
(Section  2.2).  The  details  of  the  theory  and  original  experimental  veri¬ 
fication  are  given  in  Reference  136.  This  method  has  subsequently  been  exten¬ 
sively  tested  against  newer  experimental  data  end  found  to  be  accurate  by 

Peterson1  ,  Sorrells,  Jackson  and  Czarnecki1  ,  Stallings  and  Lamb1  ,  and 
141 

Tendeland  for  a  wide  variety  of  geometries,  flow  conditions,  and  heat 
transfer  rates.  The  method  has  been  programmed  for  many  computers  and  forms 
the  skin  friction  module  of  the  Integrated  Supersonic  Design  and  Analysis 

Q7 

System' ' . 

Consistent  use  of  the  Sommer  and  Short  method  from  the  preliminary  design 
phase  through  flight  testing  will  improve  drag  correlations.  Note  that  drag 
component  correlations  used  during  the  preliminary  design  phase  are, 
ultimately,  based  on  older  wind  tunnel  and  flight  test  data.  A  considerable 
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amount  of  these  data  were  reduced  using  the  Sommer  and  Short  method  as  a 
basis.  It  is  consistent  to  continue  use  of  this  method  for  subsequent  drag 
buildups  on  new  configurations. 

4.2  Wave  Drag 

The  wave  drag  produced  by  standing  pressure  waves  that  are  not  possible 

5 

in  subsonic  flow  has  been  traditionally  decomposed  into  two  components  :  a 

zero-lift  wave  drag  and  a  wave  drag-due-to-1 ift.  The  decomposition  is  a 

result  of  the  supersonic  area  rule  method  used  for  wave  drag  prediction  and  is 

not  permitted  if  nonlinear  methods  are  used.  These  drag  components  vary  as: 

Zero  Lift:  D  (volume)2/(length)^ 

o  wave 

Lift-dependent  or  induced:  waye  (M2-l)(lift)2/(lifting  length)2 

and  must  be  added  to  the  supersonic  counterparts  of  the  subsonic  skin  friction 
drag  and  induced  or  vortex  drag.  Thus,  the  supersonic  design  problem  of 
obtaining  maximum  lift-to-drag  ratio  is  more  complicated  due  to  constraints 
arising  from  these  additional  drag  terms. 

Linear  supersonic  aerodynamics  methods  are  the  mainstay  of  the  aircraft 
industry  and  are  routinely  used  for  preliminary  design  because  of  their 
simplicity  and  versatility  in  spite  of  their  limitations  to  slender  configura¬ 
tions  at  low  lift  coefficients.  Not  surprisingly,  most  successful  supersonic 
designs  to  date  have  adhered  to  the  theoretical  and  geometrical  limitations  of 
these  analysis  methods.  Linearized  theory  has  not  been  supplanted  since 
second-order  theories  have  not  been  as  successful  at  predicting  experimental 
data,  ar.d  practical  application  of  nonlinear  methods  has  been  precluded  by 
severe  geometrical  restrictions,  i.e.  axisymmetric  bodies,  conical  flow,  etc. 

The  normal  pressure  drag  components,  consisting  of  form  drag,  vortex 
drag,  lift-dependent  wave  drag  and  volume-dependent  wave  drag,  are  shown  in 
Fig.  b. 
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The  supersonic  form  d  r  ci  g ,  or  drag  due  to  adding  the  boundary-layer 

displacement  thickness  to  the  physical  areas  to  obtain  the  effect  of  the 

boundary  layer  on  the  wave  drag,  is  conventionally  ignored^’ 

Calculations  have  shown  this  contribution  is  small  (±1%)  and  dependent  upon 

19 

longitudinal  area  distribution  for  bodies  of  revolution  .  This  assumption  is 

55 

well  within  the  accuracy  of  the  prediction  capability  for  the  body  alone  * 

^ ,  particularly  for  low  fineness  ratio  bodies^.  The  accuracy  of  this 
assumption  is  unknown  for  wings. 

As  shown  in  Fig.  28,  both  near-  and  far-field  methods  are  used  to 

calculate  the  components  of  inviscid  pressure  drag,  Fig.  6. 

39 

A  surface  panel  method,  such  as  PAN  AIR  ,  that  represents  the  config¬ 
uration's  thickness,  camber  and  angle  of  attack  with  surface  distributions  of 
sources  and  doublets  and  sums  the  resulting  pressure  times  the  streamwise 
component  of  area  (near-field  method)  obtains  the  sum  of  the  vortex  drag, 
lift-dependent  and  volume-dependent  wave  drags.  The  use  of  these  methods 
offer  improved  geometric  representation  of  the  configuration,  but  do  not 
permit  independent  consideration  of  each  drag  component  for  optimization 
purposes. 

Far-fielrf  methods  relate  the  wave  drag  to  the  lateral  convection  of 

streamwise  momentum  through  a  cylindrical  control  surface,  Fia.  28,  placed  a 

large  distance  (several  wingspans)  from  the  conf iyuration.  The  vortex  drag  is 

similarly  related  to  the  transverse  components  of  momentum  through  a  crossflow 

plane  far  downstream  in  the  wake  --  the  Trefftz  plane.  The  supersonic  vortex 

drag  is  identical  to  the  induced  drag  for  subsonic  flow  since  the  trailing 

vortic.ity  remains  essentially  stationary  with  the  fluid,  regardless  of  flight 

speed.  Due  to  the  inherent  theoretical  limitations  of  each  approach,  a 

combination  of  near-  and  far-field  methods,  along  with  semiempirical  modifica- 

5  97 

tions,  are  currently  used  for  drag  prediction  and  design  optimization  * 
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The  zero-lift,  or  volume-dependent,  wave  drag  is  commonly  computed  using 
a  numerical  implementation  of  the  Supersonic  Area  Rule  concept,  proposed 
somewhat  independently  by  Hayes,  Jones,  and  Whitcomb.  This  far-field  method 
relates  the  zero-lift  wave  drag  of  a  wing,  fuselage,  empennage,  nacelles, 
etc.,  to  a  number  of  developments  of  the  normal  components  of  cross-sectional 
areas  as  intersected  by  inclined  Mach  planes  through  the  Von  K&rm&n  slender- 
body  formula  (Fig.  28).  The  use  of  linear  source-sink  distributions  to 
represent  the  configuration  precludes  representation  of  both  local  lift  and 
total  integrated  lift.  The  wave  drag  due  to  lift  and  the  vortex  drag  are, 
perforce,  calculated  by  another  method,  conventionally  the  near-field  Mach  Box 
method  for  thin  (zero  thickness)  wings. 

With  an  additional  restriction  on  the  slenderness  of  the  configuration, 
the  Transfer  Area  Rule  can  be  obtained  that  permits  optimization  of  the 
wing-body  for  minimum  wave  drag  due  to  volume,  independent  of  the  wing.  Since 
the  fuselages  of  military  aircraft  are  usually  not  slender,  the  use  of  the 
Transfer  Area  Rule  in  this  case  would  represent  a  greater  violation  of  theory 
than  use  of  the  Supersonic  Area  Rule. 

Numerous  comparisons  of  this  theory  with  experimental  data  from  the 

5  10 

slender  XB-70  and  Supersonic  Transport  models  ’  have  demonstrated  that  it 
is  remarkably  accurate  (less  than  ±  10%  error  in  total  drag)  over  a  range  of 
Mach  numbers  from  1.2  to  3.0. 

Bonner^0,  using  an  integral  similar  to  Von  Kirm&n's  slender  body  formula 
derived  by  Hayes,  obtains  both  the  zero-lift  and  lift-dependent  wave  drag. 

This  procedure,  however,  requires  a  priori  knowledge  of  the  longitudinal 
distribution  of  lift.  The  use  of  this  procedure  is  relatively  limited  since  a 
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near-field  method  capable  of  obtaining  both  the  vortex  and  lift-dependent  wave 
drags  must  be  used  to  supply  the  required  lift  distribution  and  lifting 
surfaces  are  much  less  likely  to  satisfy  the  slenderness  requirements .  The 
distant  point  of  view  does,  however,  incorporate  lift-volume  interference  not 
reflected  in  many  of  the  theoretical  techniques  used  to  estimate  the  surface 
pressure. 

4.3  Lift-Induced  Drag 

The  drag  due  to  lift,  both  wave  and  vortex,  has  conventionally  been 
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calculated  using  the  Mach  Box  method  of  Middleton  and  Carlson'  ,  currently 
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with  many  semi  empirical  modifications  to  alleviate  known  shortcomings  .  A 
comprehensive  review  of  this  method,  as  of  1974,  is  contained  in  Ref.  23  and 
later  modifications  are  contained  in  Ref.  21  ana  Ref.  97.  The  method  sub¬ 
divides  the  zero-thickness  planform,  often  with  the  fuselage  outline  (Fig. 

28),  into  rectangular  panels  whose  diagonals  are  aligned  with  the  free-stream 
Mach  lines  of  the  flow.  Alternatively,  the  sides  of  the  element  could  be 
parallel  to  the  Mach  lines  and  the  Mach  Diamond  or  Characteristic  Box  method 
would  result.  Generally,  the  Mach  Box  method  is  more  accurate  for  supersonic 
edges,  while  the  Characteristic  Box  gives  better  results  for  wings  with 
subsonic  edges®1. 

The  linear  theory  integral  equation  relating  camber  surface  slope  to 
lifting  pressure  difference,  AC^,  across  the  planform  is  numerically  evaluated 
by  considering  ACp  a  constant  within  each  box  or  element  and  reducing  the 
integral  equation  to  an  algebraic  summation.  The  summation  can  be  solved  for 
either  the  camber  surface  slope  for  a  known  pressure  distribution  (wing  design 
case)  or  the  lifting  pressure  coefficient  in  terms  of  known  camber  surface 
slopes  (lift  analysis  case).  A  number  of  component  loadings  (up  to  seventeen) 
can  be  combined  in  the  design  mode  to  determine  optimum  camber  shapes  for  minimum 
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drag  due  to  lift  at  a  given  lift  coefficient,  subject  to  optional  pitching 

moment  constraints.  The  effects  of  the  fuselage  and  nacelles  may  be  included 

in  the  component  loadings  and  additional  constraints  may  be  applied  to  the 

design  pressure  distribution  and  local  camber  surface  to  provide  physical 

realism.  Linear  theory  consistently  overestimates  expansion  pressures  (to 

1 14 

less  than  vacuum  pressure!)  and  underestimates  the  compression  pressure 

Proper  selection  of  the  component  loadings  and  constraints  requires  complete 

understanding  of  the  computer  programs  backed  by  critical  operational  use, 
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otherwise,  occasional  "endless  loops"  may  result  .  There  may  also  be  small 
discrepancies  between  wing  loadings  and  forces  determined  for  an  optimized 
wing  and  the  leadings  and  forces  for  the  same  shape  upon  submittal  to  the 
evaluation  program  . 

The  preceding  supersonic  aerodynamic  design  and  analysis  methods  have 
been  substantiated  for  numerous  wings,  wing/body  and  wing/body/nacelle  config¬ 
urations  at  the  lower  lift  coefficients  typical  of  supersonic  cruise  condi¬ 
tions.  A  comparison  with  wind-tunnel  data  from  the  Lightweight  Experimental 
Supercruise  Model  (LES)  at  the  design  Mach  number  and  lower  lift  coefficients 
is  shown  in  Fig,  30.  Test-theory  comparisons  for  supersonic  transport  and 
bomber  models  at  typically  higher  cruise  Mach  numbers  and  lower  design  lift 
coefficients  are  equally  as  accurate^’  although  only  if  performed  by 
experienced  designers  familiar  with  the  methods.  Fig.  31. 

The  major  unresolved  drag  prediction  problem  of  thin,  highly  swept  wings 
typical  of  efficient  supersonic  flight  is  the  evaluation  of  leading-edge 
thrust  or  leading-edge  suction.  The  large  influence  on  the  drag  polar  of  the 
various  assumptions  used  to  determine  the  magnitude  of  this  force  is  shown  in 
Fig.  32  for  the  same  configuration  as  shown  in  Fig.  30.  The  polar  shown 
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in  Fig.  32  is  for  M  =  1.35,  an  off-design  condition  for  the  M  =  1.8  super¬ 
cruiser  configuration. 

The  leading  edge  thrust  results  from  the  low  pressure  induced  by  the  high 
flow  velocities  around  the  leading  edge  from  the  stagnation  point  on  the 
undersurface  of  the  wing  to  the  upper  surface.  For  the  high  aspect  ratio 
wings  typical  of  low  subsonic  flight  speeds,  this  force  largely  counteracts 

OO  n  t  c 

the  drag  from  the  pressure  forces  acting  on  the  remainder  of  the  airfoil4-'"’ 

The  thrust  force  diminishes  with  increasing  speeds,  but  exists  as  long  as  the 

leading-edges  are  subsonic  .  It  was  found  to  be  negligible  at  M  =  3.0  for 

the  configuration  analyzed  by  Kulfan^.  Thus,  thrust  effects  have 

conventionally  been  ignored  for  supersonic  wing  design  and  analysis  at  Mach 

numbers  approaching  three^6,  with  excellent  results,  as  shown  in  Fig.  31. 

The  renewed  interest  in  this  phenomena  resulted  from  the  increased 

sophistication  in  supersonic  wing  design  required  to  achieve  higher  design 

lift  coefficients  for  the  lower  Mach  numbers  typical  of  proposed  fighter 
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designs  with  sustained  supersonic  cruise  capability  .  The  maneuver  require¬ 
ments  also  demand  higher  available  lift  coefficients  and  the  attendant  thick¬ 
ness  and  complex  camber  shapes  to  minimize  drag. 

The  current  complex  of  NASA  computer  programs^'*  contains  four  user 
selected  options  to  estimate  leading-edge  thrust,  Fig.  32.  The  basic  Mach  Box 
zero-thickness  wing  design  and  analysis  method  provides  the  zero  leading-edge 
suction  drag  estimate  since  the  low  pressures  do  not  have  a  forward  facing 
area  to  act  upon.  The  Polhamus  analogy  assumes  the  flow  is  separated  from  the 

entire  leading-edge  and  has  formed  the  spiral  vortex  sheets  characteristic  of 

76 

slender  wings  at  higher  angles  of  attack  .  The  Polhamus  analogy  has  been 

7  c 

extended,  by  Kulfan  ,  to  account  for  airfoil  shape,  wing  warp  and 
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plan-form  effects  in  order  to  determine  the  influence  of  wing  geometry  on 
leading-edge  vortices.  He  demonstrates  that,  except  for  very  slender  wings, 
lowest  drag  due  to  lift  is  achieved  with  attached  flow. 

The  attainable  thrust  option  (Fig.  32),  developed  by  Carlson,  Mack  and 
Barger  ,  has  been  found  to  work  well  for  wings  with  standard  NASA  airfoil 
sections,  from  which  the  correlations  were  developed1  .  However,  the  attain¬ 
able  thrust  forces  predicted  did  not  agree  with  those  experimentally  measured 

on  wings  with  sharp  arid  varying  leading-edge  radii;  large  differences  in  both 
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levels  and  trends  were  evident 

The  magnitude  of  the  full  leading-edge  thrust  is  dependent  upon  the 
upwash  in  the  vicinity  of  the  leading  edge  and  the  effect  on  the  pressure 
distribution.  It  is  determined  by  calculating  the  limit  of  the  pressure 
coefficient  (which  theoretically  approaches  infinity)  multiplied  by  the 
distance  from  the  leading  edge.  The  evaluation  of  this  limit  is  difficult 
because  the  pressure  distributions,  as  calculated  by  panel  methods,  vary  from 
the  exact  theoretical  values  required.  This  limit  problem  and  the  resulting 
impact  on  wing  design  and  drag  prediction  have  been  extensively  studied  by 

pi 

Carlson  and  Mack^'  . 

Calculations  using  near-field  surface-panel  methods  should  predict  the 
full  suction  drag  polar;  however,  the  pressure-area  integral  is  sensitive  to 
leading-edge  paneling  density.  Near  thin  wing  tips,  the  paneling  may  be 
insufficient  for  the  accurate  evaluation  of  this  integral  and  the  resulting 
force  will  be  bounded  between  the  full  and  zero  leading-edge  suction  values. 
The  force  will  thus  depend  ori  the  numerical  representation  of  the  wing,  rather 
than  aerodynamic  assumptions,  a  disconcerting  situation  for  the  designer. 
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As  observed  by  several  investigators,  a  reasonable  estimate  for  the  drag 
polar  could  be  obtained  by  numerically  averaging  the  zero  and  full  leading- 
edge  suction  drag  polars. 

Drag  polars  calculated  by  two  near-field  methods  for  Mach  numbers  of 
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1.196,  1.6,  and  2.2  are  shown  in  Fig.  33  .  The  complex  canard-wing-fuselage 

configuration  and  paneling  arrangement  are  shown  in  Fig.  21.  The  skin 

friction  drag  was  estimated  using  the  methods  of  Section  4.1.  Considerable 

scatter  is  evident  in  the  zero  lift  drag  which  may  be  attributed  to  either  the 

skin  friction  estimate  or,  more  likely,  the  zero  lift  wave  drag  predictions  of 

the  near-field  methods.  The  variation  in  vortex  and  lift-dependent  wave  drag 

with  angle  of  attack  was  correctly  predicted  by  the  higher-order  panel  method, 

PAN  AIR^.  The  lower-order  method,  USSAERO-B,  has  recently  been  updated  to 

improve  supersonic  prediction  capability  by  inclusion  of  the  Triplet 
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Singularity  *  .  The  wing  camber  design  Mach  number  for  this  configuration 

was  1.8.  Many  additional  test-theory  comparisons,  exhibiting  much  closer 
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agreement,  are  contained  in  Tinoco,  Johnson,  and  Freeman  and  Tinoco  and 
146 

Rubber!  .  The  prediction  methods  included  a  pilot  version  of  PAN  AIR,  the 
Mach  Box  method  and  the  Woodward  Constant-Pressure  panel  method  (USSAERO)  for 
a  Mach  three  Recce-Strike  design  and  the  supercruise  configuration  shown  in 
Figs.  30  and  32. 

Correction  procedures  to  remedy  the  observed  shortcomings  of  linear 
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theory  are  under  continuous  development.  Stancil  has  modified  near-field 
linear  theory  to  improve  the  prediction  of  compression  and  expansion  pressures 
and  to  eliminate  local  singularities.  His  modifications,  the  use  of  exact 
tangency  surface  boundary  conditions  and  the  corrected  local  value  of  Mach 
number,  cause  the  solution  technique  to  become  iterative.  The  modifications 
improve  the  zero-lift  wing  wave  drag  prediction  capability  of  linear  theory 
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for  the  F-8C  .  Shrout  and  Covel!  compared  both  Stancil's  modified  linear 

theory  and  the  far-field  slender-body  theory  zero  lift  wave  drag  predictions 

with  their  experimental  data  for  a  series  of  forebody  models  having  various 

degrees  of  nose  droop.  The  far-field  method  correctly  predicted  trends  for 

Mach  numbers  of  1.2  and  1.47,  but  was  erroneous  at  1.8  for  the  cambered 
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models.  Stancil's  method  '  correctly  predicted  wave  drag  trends  for  Mach 
numbers  of  1.47,  1.80  and  2.16. 

Another  breakdown  of  linear  theory  occurs  when  the  crossflow  velocities 
exceed  sonic  value,  a  situation  encountered  on  fighter  aircraft  at  maneuver 
lift  coefficients.  The  formation  of  crossflow  shock  waves  causes  the  experi¬ 
mental  drag-due-to-1 i ft  performance  of  wings  designed  by  linear  theory  to  fail 
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to  meet  predicted  values.  The  supercritical  conical  camber  concept  ,  (SC  ), 
has  been  proposed  to  regain  wing  performance  for  wings  swept  such  that  the 
Mach  number  normal  to  the  leading  edge  is  transonic.  The  design  method 
employs  repetitive  application  of  a  transonic  full  potential  flow  solver, 
COREL,  and  a  specially  adapted  version  of  the  Woodward  USSAERO  linear  theory 
paneling  code  .  Comparison  of  this  design  procedure  with  data  from  the  wind 
tunnel  tests  of  a  demonstration  wing  show  that  substantial  reductions  in 

Q 1 

drag-due-to-lift  are  possible,  Fig.  34,  at  the  maneuver  design  condition  . 

Bonner  and  Gingrich1^  approach  fighter  wing  design  through  the  use  of 
variable  camber  and  a  design  compromise  between  the  transonic  maneuver  config¬ 
uration  and  an  optimum  supersonic  cruise  design.  The  unconstrained  supersonic 
design  was  initially  accomplished  by  conventional  linear  theory  while  both 
small -disturbance  and  full -potential  theories  were  used  for  the  transonic 
design.  These  designs  established  the  boundaries  for  the  design  cycle. 
Experience  then  guided  the  selection  of  the  compromise  camber  shape,  followed 
by  alternating  constrained  optimizations  at  the  supersonic  and  transonic 
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design  points.  This  multipoint  design  process  illustrates  the  need  for 
accurate  drag  prediction.  Spurious  high  drag  regions  entirely  due  to  the  use 
of  linear  theory,  shown  in  Fig.  35  at  the  sonic,  leading-edge  condition,  could 
eliminate  the  optimum  design  from  consideration  on  the  first  cycle! 

Significant  potential  for  further  improvement  in  wing  design  for  minimum 
drag-due-to-lift  exists  through  the  use  of  nonlinear  computational  methods 
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early  in  the  design  cycle  .  A  comprehensive  set  of  design  constraints  for 

7 

efficient,  highly  swept  wings  has  been  determined  by  Kulfan  . 

The  problems  of  fighter  design  and  analysis  are  further  compounded  by  the 

fact  that  linear  theory  does  not  apply  to  relatively  large  fuselages  with 
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nonsmooth  area  distributions,  as  shown  by  Wood,  et.al.  Four  wings  were 

designed  and  tested  with  a  common  fuselage  in  an  effort  to  adapt  supersonic 

technologies  to  a  M  =  1.8  fighter  design.  Their  test-theory  comparisons,  Fig. 

36,  clearly  show  the  inaccuracy  in  zero-lift  drag  prediction  using  the 
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Middleton  and  Lundy  complex  of  computer  programs.  This  fuselage  severely 
violates  the  assumptions  of  linear  theory  concerning  smooth,  slender  area 
distribution  bodies.  Further  inaccuracies  as  the  configuration  is  built  up 
are  opposite  in  sign  and  cancel  most  of  the  fuselage  error  —  illustrating  the 
somewhat  forgiving  nature  of  linear  aerodynamics.  The  zero-lift  drag  compari¬ 
son  for  each  of  the  four  wings  is  shown  in  Fig.  37.  The  prediction  capa¬ 
bilities  of  linear  theory  are  excellent  for  the  complete  configurations  and 
vary  with  slenderness  and  smoothness  of  the  area  distribution.  The  drag-due- 
to-lift  comparisons  were  equally  as  accurate;  however,  they  were  performed 
using  in-house  empirical  corrections. 
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Linear  theory,  as  developed  by  Middleton  and  Lundy  ,  has  become  the 
standard  preliminary  design  tool  for  supersonic  aircraft.  Correction  proce¬ 
dures  to  alleviate  the  well-known  restrictions  of  the  theory  are  under 
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constant  development.  Higher-order  panel  methods,  necessary  for  supersonic 
flow  analyses,  along  with  nonlinear  numerical  methods,  are  also  used  during 
the  latter  design  stages,  particularly  for  lower  design  Mach  number,  blunt 
fighter  configurations. 

An  innovative,  statistical  method  for  the  prediction  of  fighter  aircraft 
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drag-due-to-lift  was  developed  by  Simon,  et.al.  ,  for  Mach  numbers  of  0.4 

through  2.5.  A  large  experimental  data  base  was  fitted  with  regression 

equations  based  on  easily  calculated  nondimensional  geometric  parameters.  The 

accuracy  of  the  resulting  regression  equations  was  only  limited  by  the  expanse 
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of  the  data  base  from  which  they  were  obtained,  as  shown  by  Johnson  and 
Tomasetti 

5  Numerical  Aerodynamics 

Here,  the  term  numerical  aerodynamics  is  used  to  delineate  the  full 

spectrum  of  aerodynamic  methods  requiring  large  computers  —  from  the 

linearized  theory  panel  methods  to  the  full  Navier-Stokes  equations.  The  term 

computational  fluid  dynamics  (CFD)  is  often  used  in  this  context,  and  just  as 

often  limited  to  mean  only  Navier-Stokes  computations,  particularly  with 

respect  to  future  goals.  Numerical  aerodynamics  now  permeates  nearly  every 

stage  of  the  design  process,  as  well  as  each  speed  regime. 

Many  prognostications  concerning  the  possible  roles  of,  primarily,  CFD 

and  the  wind  tunnel  have  been  published  (Ref.  27,  29,  30,  42,  48,  53,  73,  79, 

83,  99,  110,  115,  121,  124,  125,  132,  167)  since  Chapman's  1975  article2®  and 

the  prediction  that  CFD  would  supplant  traditional  wind  tunnel  tests  within  a 

decade.  The  pacing  item  in  CFD  technology  at  that  time  was  computer  size  and 

30 

speed.  A  later  article  recognized  turbulence  models  to  be  a  pacing  item. 

In  addition  to  the  long-term  problem  of  turbulence  modeling,  the  problem  of 
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proper  grid  generation  is  of  immediate  concern  to  those  engaged  in  the  use  of 

CFD  codes  in  the  production  (as  apposed  to  research)  mode.  Computer  hardware 
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has  advanced  tremendously  in  the  last  decade  ,  as  have  flow  simulations  using 

the  Reynolds  averaged  Navier-Stokes  equations.  Still,  the  wind  tunnel  is 

the  primary  tool  for  design  verification  and  a  large  portion  of  design 

development  work.  CFD  analysis  of  full  three-dimensional  configurations 

capable  of  flight  may  still  be  a  decade  away^’  The  fundamental 

limitation  of  CFD  during  the  last  decade  --  turbulence  modeling  --  remains  the 

pacing  item.  Progress  in  this  area  in  the  next  15  years  is  not  projected  to 

remove  this  limitation'  and  recent  reviews  of  turbulence  models  9  *  5 
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9  reinforce  this  contention.  The  fundamental  physics  of  turbulence  is 
not  known  well  enough  to  permit  accurate  mathematical  modeling  and  flow 
simulation.  If  the  present  is  any  indication  of  the  future,  drag  prediction 
and  design  optimization  with  drag  as  an  objective  function  using  CFD  are  even 
further  in  the  future. 

Recent  publications  recognize  the  complementary  roles  of  wind  tunnels  and 
CFD  ,  while  the  role  of  the  lower-order  approximations  (linear,  potential, 
Euler)  is  often  understated.  Each  of  these  methods  will  be  found  in  a  comple¬ 
mentary  place  in  the  design  process  depending  on  configuration,  schedule, 
cost,  availability  of  hardware  and  software,  and  the  persuasion  of  the 
designer  and  his  management.  Confidence  in  the  accuracy  of  the  predictions  is 
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tantamount  .  Aircraft  design  departments  should  evolve  a  full  spectrum 

of  numerical  aerodynamics  methods,  coupled  through  a  data  management  system, 
to  provide  accurate,  timely,  cost-effective  design  data1  9  L  9  9  "  °. 
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CLEAN  WING  REFERENCE 
CONFIGURATION 
•  Isolated  wing/body  modal 


AERODYNAMIC  REFERENCE 
CONFIGURATION 

•  Wlng/body  model  with  flow-throunh 
nacallai  and  reference  nonles 


PROPULSION  REFERENCE 
CONFIGURATION 

•  Wing/body  model  with  blown 
nacelles  end  reference  nonles 


Figure  1.  Test  Approach 
(from  Tjonneland,  et  al.  Ref.  147) 
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Figure  2.  Transport  Aircraft  Drag  Buildup 
(Full  Scale,  M  =  0.8) 


3.0 


igure  4.  Drag  Data  Scatter  About  Drag  Polar  Fairings  with  Mach  Numbe 

(From  Arnaiz,  Ref.  3) 
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Figure  5.  Typical  Aircraft  Drag  Breakdowns 
(from  Butler,  Ref.  19) 
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Figure  6.  Components  of  Aircraft  Drag 
(from  Butler,  Ref.  19;  See  also  Leyman,  Ref.  81) 
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Figure  8.  Subsonic  Parasite  Drag-Bomber /Transport 
(frdm  Paterson,  et  al,  Ref.  109) 


Figure  9.  Subsonic  Parasite  Drag-Fighter/Attack 
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FigurG  12,  Equivalent  Body  Concept 
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Figure  14. 
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of  Empirical  Flat  Plate  Skin  Friction 
for  Incompressible  Turbulent  Flow. 
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(from  MacWilkinson,  et  al 
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Figure  19.  Lift  and  Speed  Regions  for  Calculation  of  Drag  Due  to  Lift 

(from  Schemensky,  Ref.  122) 
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is  the  minimum  profile  drag  =  Cp0  -  KClo 

is  the  zero  lift,  drag  coefficient  _ 


is  the  theoretical  drag  due  to  lift  = 
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■  CT  )  is  the  additional  profile  drag 
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is  the  lift  coefficient  for  drag  break,  also  f-Lpu 

is  the  intermediate  drag  break  increment,  also  Z\Cpp 

is  the  drag  variation  for  C  >CT 
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Figure  20.  The  Subsonic  Drag  Polar 
(from  Simon,  et  al.  Ref.  131) 


COEFFICIENT  OF  LIFT 


(from  Chen,  et  al.  Ref.  31) 


Figure  23.  Historic  Correlation  of  Transonic  Drag  Rise 
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CD 


Figure  25.  F-5  versus  C^;  M  =  0.9 

(from  Axelson,  Ref.  4) 
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Figure  27.  Correlation  of  Calculated  and  Ileasured  Wing-Body 

Drag  Divergence  ?lach  Number,  Compressibility  Drag, 
and  Drag  Rise  for  Three  Supercritical  Wings. 

(from  Henne,  et  al.  Ref.  57) 
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Figure  29.  Supersonic  Skin  Friction 
(from  Leyman,  et  al.  Ref.  81) 
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TEST-THEORY  COMPARISON 
LES-216  WIND  TUNNEL  MODEL 


Figure  30.  Supercruiser  Type  Configuration  Test-Theory  Comparison 

(from  Kulfan,  et  al,  Ref.  78) 
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igure  31.  Experimental  and 
Supersonic-Cru 

(from  Robins,  et 


CL 


[  Theoretical  Drag  Polars  of 
tise  Airplanes,  Rc  =  4.8  x  106 


al,  Ref.  116) 


Figure  32.  Leading  Edge  Suction  Effects,  LES  216 
(from  Middleton,  et  al,  Ref.  97) 
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Figure  33.  Drag  Polars,  CDAF  Configuration 
(from  Miller,  et  al.  Ref.  98) 
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Figure  35.  Supersonic  Volumetric  Efficiency 
(from  Bonner,  et  al.  Ref.  12) 


NOTE:  PERCENTAGES  ARE  REFERENCED  TO 
WIND  TUNNEL  DATA 


Figure  37.  Zero  Lift  Drag  Comparison 
(from  Wood,  et  al.  Ref.  164) 


